This is the second part of the analysis. In the first part (see
../input/PHO5-data/20231019-pool-qc-PHO5.Rmd), I did QC and
exported the filtered dataset. Here, I will continue working with that
dataset and answer our biological questions.
Data
Import the background subtracted data
dat0 <- read_tsv("../input/20231023-PHO5-bg-subtracted-data.tsv", col_types = "ccccdddddc")
Filter the data
dat <- filter(dat0, host != "PHO84", flag == "pass", date != "02/10") %>%
# based on previous QC, the following sample (both replicates) have high
# variance - one biological replicate is highly expressed, while the other
# two have mNeon, but barely any RFP expression.
mutate(
host = fct_recode(host, pho2 = "pho2∆"),
flag = ifelse(plasmid == "233" & host == "pho2", "high.var", flag))
Number of replicates left for each sample
expt <- dat %>%
filter(host %in% c("PHO2", "pho2"), !plasmid %in% c("188", "194")) %>%
group_by(date, plasmid, host) %>%
summarize(n = n(), .groups = "drop")
expt %>%
ggplot(aes(x = plasmid, y = n)) +
geom_col(aes(fill = host)) +
facet_grid(date ~ .) +
scale_fill_manual(values = c("PHO2" = "gray30", "pho2" = "gray70")) +
theme_minimal() + background_grid(major = "none") + panel_border(size = 0.5) +
scale_y_continuous(name = "Replicates", breaks = c(6)) + xlab(NULL) +
theme(axis.text.x = element_text(angle = 90),
strip.text.y = element_text(angle = 0),
legend.position = "top")

Chimera makeup information
meta <- read_tsv("../input/20230208-chimera-Pho4-makeup.txt", col_types = "ccccc")
Summarize data
Here we would like calculate the ratio of RFP/GFP for each chimera
(plasmid) across all replicates, including from different days. Note
that the parameter of interest is a ratio, which can be estimated using
either “means of ratios” or “ratios of means”. These are just two
specific instances of a more general estimator, representing two choices
of the weights. The “means of ratios” first calculates the ratios for
each replicate within a plasmid, then average them. In this calculation,
each replicate is given the weight of 1/n (equal). The “ratios of means”
first sum up the GFP and RFP values separately across the replicates for
each plasmid, then take the ratio between them. In this estimator, the
weight for each replicate is x / sum(x), where x is the denominator in
the ratio, i.e., GFP. In other words, this estimator will give more
weights to the replicates where the chimera had a higher expression
level.
Both estimators are known to be biased. We will ignore that for the
moment. In terms of a choice between the two, it seems that there is no
reason to give more weights to the experiments with a higher GFP signal.
So, the “means of ratios” seems a more natural choice. However, we will
calcultae both and decide later.
A final question is how to calculate the variance of the ratio
estimate. According to the survey package manual,
an approximate estimator for the variance is
\[
r = \frac{\bar{y}}{\bar{x}}, \text{where}\
\bar{y}=\frac{1}{n}\sum_{i=1}^{n}y_i\ \text{and}\
\bar{x}=\frac{1}{n}\sum_{i=1}^{n}x_i\ \\
\hat{V}(r) = (1-\frac{n}{N})(\frac{1}{\bar{x}^2})\frac{s_r^2}{n}\
\text{where}\ s_r^2=\frac{1}{n-1}\sum_{i=1}^{n}(y_i-rx_i)^2
\]
Assuming that N>>n, we can ignore the first term in the
variance estimator. The rest can be calculated from the data
datsum <- dat %>%
filter(!is.na(plasmid)) %>%
group_by(plasmid, host) %>%
summarize(
n = n(),
mG = mean(BL1.H),
mR = mean(YL2.H),
A = mean(YL2.H/BL1.H),
r = mR/mG,
s2 = 1/(n-1)*sum((YL2.H - r*BL1.H)^2),
vr = 1/(mG^2)*s2/n,
se = sqrt(vr),
.groups = "drop"
) %>%
select(-s2, -r, -vr)# %>%
#pivot_wider(names_from = host, values_from = BL1.H:`nR/G`) %>%
#mutate(`pho2∆/PHO2` = `R/G_pho2∆`/`R/G_PHO2`,
# `n.pho2∆/PHO2` = `nR/G_pho2∆`/`nR/G_PHO2`)
For each chimera, we would also like to calculate three
values:
- A in pho2∆: this is its base activity without Pho2
- A in PHO2: this is its full activity with Pho2
- A_PHO2 / A_pho2∆: this is the Pho2 enhancement of activity
We assign the chimeras into several groups, based on their A_PHO2 and
A_PHO2/A_pho2∆
ximera <- datsum %>%
pivot_wider(id_cols = plasmid, names_from = host,
values_from = c(A, se)) %>%
mutate(
rA_PHO2 = A_PHO2 / A_PHO2[plasmid == "194"],
rA_pho2 = A_pho2 / A_pho2[plasmid == "194"],
boost = A_PHO2 / A_pho2,
group = case_when(
plasmid %in% c("188", "194") ~ "ref",
rA_PHO2 < 0.2 ~ "n.f.",
.default = "chimera"
),
group = fct_relevel(group, "ref", "chimera", "n.f.")
) %>%
left_join(select(meta, plasmid, set, symbol), by = "plasmid") %>%
mutate(symbol = fct_reorder(symbol, rA_PHO2, .desc = TRUE)) %>%
relocate(c(set, symbol, group), .after = plasmid)
Export the summarized data
write_tsv(ximera, file = "../output/20231125-PHO5pr-chimera-summarized.tsv")
To be able to plot all the data points, let’s generate another data
frame with the individual ratios.
dat_sep <- dat %>%
filter(!is.na(plasmid)) %>%
mutate(A = YL2.H/BL1.H) %>%
select(plasmid, host, BL1.H, YL2.H, A, flag) %>%
left_join(select(meta, plasmid, set, symbol), by = "plasmid")# %>%
#mutate(symbol = fct_reorder(symbol, rA_PHO2, .desc = TRUE))# %>%
Analysis
Plotting functions
—>
Source the scripts
my_plot_ratio <- function(selection){
# custom colors for this function
date.colors = c(brewer.pal(name="Dark2", n = 8), brewer.pal(name="Paired", n = 8))
host.colors = c("PHO2" = "gray30", "pho2" = "gray70")
point.colors = c("PHO2" = "forestgreen", "pho2" = "purple4")
# prepare data
tmp <- my_data_prep(selection)
# plotting
p <- tmp %>%
select(-c(FSC.H, nGFP, nRFP, flag)) %>%
mutate(`R/G` = YL2.H/BL1.H) %>%
pivot_longer(cols = c(BL1.H, YL2.H, `R/G`),
names_to = "parameter", values_to = "value") %>%
mutate(parameter = factor(parameter, levels = c("R/G", "YL2.H", "BL1.H"),
labels = c("RFP/GFP", "PHO5pRFP", "Pho4-GFP"))) %>%
ggplot(aes(x = symbol, y = value, group = host)) +
stat_summary(aes(group = host), fun.data = "mean_cl_boot", geom = "errorbar",
position = position_dodge(0.5), width = 0.3) +
geom_bar(aes(fill = host), width = 0.5, alpha = 0.8,
stat = "summary", fun = "mean", position = position_dodge(0.5)) +
geom_point(data = function(x) subset(x, !symbol %in% c("CCCCC", "SSSSS")),
aes(group = host, color = date), size = 1, shape = 3, alpha = 0.9,
position = position_jitterdodge(dodge.width = 0.5, jitter.width = 0.1)) +
scale_color_manual(values = date.colors, guide = "none") +
#geom_point(data = function(x) subset(x, !symbol %in% c("CCCCC", "SSSSS")),
# aes(group = host, color = host), size = 1, shape = 3, alpha = 0.9,
# position = position_jitterdodge(dodge.width = 0.5, jitter.width = 0.1)) +
#scale_color_manual(values = point.colors) +
scale_fill_manual(values = host.colors) +
facet_grid(parameter~group, scales = "free", space = "free_x") +
theme_bw(base_size = 18) + background_grid(minor = "none") +
xlab("Pho4 chimera") +
theme(axis.text.x = element_text(angle = 30, hjust = 1, family = "mono"),
legend.position = "top",
axis.title = element_blank())
return(p)
}
Modify the component plotting function for special purposes
host.labels = c("PHO2", "pho2∆")
point.colors = c("PHO2" = "forestgreen", "pho2" = "purple4")
p1 <- dat %>%
filter(!is.na(plasmid)) %>%
mutate(plasmid = fct_reorder(plasmid, BL1.H, .fun = median) %>%
fct_relevel("194", "188")) %>%
ggplot(aes(x = plasmid, y = BL1.H)) +
geom_point(aes(color = host), position = position_jitter(0.1),
size = 1.1) +
scale_color_manual("Host", values = point.colors, labels = host.labels) +
scale_y_log10(breaks = c(100, 1000, 10000), expand = expansion(mult = 0.1)) +
scale_x_discrete(expand = expansion(mult = 0.03)) +
xlab("Pho4 constructs") + ylab("Pho4-mNeon (a.u.)") +
theme_cowplot() + panel_border(color = "gray30", size = 1.2) +
theme(axis.text.x = element_text(angle = 90, size = rel(0.6), vjust = 0.5),
axis.text.y = element_text(size = rel(0.8)),
axis.title = element_text(size = rel(0.9)),
axis.line = element_blank(),
legend.position = c(0.05, 0.9),
legend.direction = "horizontal",
legend.text = element_text(face = 3))
Warning: A numeric `legend.position` argument in `theme()` was deprecated in ggplot2 3.5.0.
Please use the `legend.position.inside` argument of `theme()` instead.
p1

#ggsave("../img/20240307-Pho4-chimera-protein-level-variation.png",
# width = 6, height = 3)
Pho4 chimera protein level variation
What is the distribution of Pho4 chimera protein levels? Is the
mCherry/mNeon ratio a faithful measure of the chimera’s activities?
Distribution of Pho4-mNeon levels grouped by plasmid and host.
tmp <- dat %>%
filter(plasmid == "194", host == "PHO2")
lm <- lm(YL2.H ~ BL1.H, data = tmp)
summary(lm)
Call:
lm(formula = YL2.H ~ BL1.H, data = tmp)
Residuals:
Min 1Q Median 3Q Max
-6194.6 -1651.7 -169.8 1417.1 6612.2
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 9562.8515 1114.3148 8.582 1.48e-14 ***
BL1.H 10.4668 0.7911 13.230 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2284 on 142 degrees of freedom
Multiple R-squared: 0.5521, Adjusted R-squared: 0.5489
F-statistic: 175 on 1 and 142 DF, p-value: < 2.2e-16
p2 <- tmp %>%
ggplot(aes(x = BL1.H, y = YL2.H)) +
geom_point(size = 1.5) +
stat_smooth(method = "lm", formula = y ~ x) +
xlab("Pho4-mNeon") + ylab("PHO5pr-mCherry") +
theme_cowplot() +
panel_border(color = "gray30", size = 1.2) +
theme(axis.line = element_blank(),
axis.title = element_text(size = rel(1.2))
)
p2
ggsave("../img/20240307-ScPho4-mCherry-vs-mNeon-consistent.png",
width = 4, height = 3.5)

All chimera (ScPho4 and CgPho4 removed) protein levels by host
host.colors = c("PHO2" = "gray60", "pho2" = "orange")
tmp <- dat %>%
filter(plasmid %in% c("188", "194"),
date %in% c("02/08", "02/11", "02/18", "02/21", "02/23", "03/31")) %>%
mutate(A = YL2.H/BL1.H,
Pho4 = factor(plasmid, levels = c("194", "188"),
labels = c("ScPho4", "CgPho4")))
p3 <- tmp %>%
ggplot(aes(x = date, y = A)) +
geom_bar(aes(fill = host), stat = "summary", fun = "mean",
position = position_dodge(0.9), alpha = 0.9) +
geom_point(aes(group = host), size = 0.6, shape = 3,
position = position_dodge(width = 0.9)) +
scale_fill_manual("Host", values = host.colors, labels = host.labels) +
scale_x_discrete(labels = 1:6) +
#stat_summary(fun.data = "mean_se", geom = "pointrange", color = "red") +
facet_grid(Pho4 ~ .) +
ylab("mCherry/mNeon") + xlab("Replicate") +
theme_cowplot() + panel_border(color = "gray30", size = 1.2)+
theme(axis.text = element_text(size = rel(0.7)),
axis.title = element_text(size = rel(1)),
axis.line = element_blank(),
strip.background = element_blank(),
legend.position = "top",
legend.title = element_text(size = rel(0.9)),
legend.text = element_text(size = rel(0.8), face = 3))
p3

#ggsave("../img/20240307-CgPho4-mCherry-vs-mNeon-consistent.png",
# width = 4, height = 3.2)
# sample size per day of experiment
tmp %>% count(date, Pho4, host) %>%
pivot_wider(names_from = host, values_from = n)
cv <- dat %>%
select(-nGFP, -nRFP) %>%
pivot_longer(FSC.H:YL2.H, names_to = "parameter", values_to = "intensity") %>%
group_by(date, plasmid, host, parameter) %>%
summarize(
n = n(),
mean = mean(intensity),
cv = sd(intensity)/mean(intensity),
.groups = "drop"
) %>%
arrange(desc(cv))
High variance samples
Summarize the background subtracted data by calculating the means and
cv for each strain.
control <- filter(dat, plasmid == "194", host == "PHO2") %>%
separate(well, into = c("row", "col"), sep = 1) %>%
droplevels()
Use the control strain (pH194 with PHO2) to identify and correct for
systematic biases
gfp.model.0 <- lm(BL1.H ~ log10(events) + date + row*col, data = control)
step(gfp.model.0)
Start: AIC=1507.62
BL1.H ~ log10(events) + date + row * col
Df Sum of Sq RSS AIC
- row:col 6 96964 3731743 1499.4
- log10(events) 1 12751 3647531 1506.1
<none> 3634779 1507.6
- date 11 3415318 7050097 1581.0
Step: AIC=1499.41
BL1.H ~ log10(events) + date + row + col
Df Sum of Sq RSS AIC
- row 3 92309 3824052 1496.9
- log10(events) 1 9919 3741662 1497.8
<none> 3731743 1499.4
- col 2 1030496 4762240 1530.5
- date 11 3409784 7141528 1570.9
Step: AIC=1496.93
BL1.H ~ log10(events) + date + col
Df Sum of Sq RSS AIC
- log10(events) 1 16207 3840259 1495.5
<none> 3824052 1496.9
- col 2 1024530 4848582 1527.1
- date 11 3423810 7247862 1567.0
Step: AIC=1495.54
BL1.H ~ date + col
Df Sum of Sq RSS AIC
<none> 3840259 1495.5
- col 2 1068825 4909084 1526.9
- date 11 3422018 7262277 1565.3
Call:
lm(formula = BL1.H ~ date + col, data = control)
Coefficients:
(Intercept) date02/09 date02/11 date02/16 date02/18
1827.3 -209.7 -259.9 -351.5 -242.3
date02/19 date02/20 date02/21 date02/22 date02/23
-319.0 -499.0 -343.8 -233.1 -439.5
date03/30 date03/31 col5 col9
-533.6 -569.1 -107.4 -211.0
gfp.model.1 <- lm(BL1.H ~ date + col, data = control)
Model for mNeon
rfp.model.0 <- lm(YL2.H ~ log10(events) + date + row*col, data = control)
step(rfp.model.0)
Start: AIC=2229.67
YL2.H ~ log10(events) + date + row * col
Df Sum of Sq RSS AIC
- row:col 6 13201086 560400308 2221.1
<none> 547199222 2229.7
- log10(events) 1 13838563 561037785 2231.3
- date 11 852447575 1399646797 2342.9
Step: AIC=2221.11
YL2.H ~ log10(events) + date + row + col
Df Sum of Sq RSS AIC
<none> 560400308 2221.1
- log10(events) 1 12713396 573113704 2222.3
- col 2 65908721 626309029 2233.1
- row 3 154001100 714401408 2250.1
- date 11 851423159 1411823467 2332.2
Call:
lm(formula = YL2.H ~ log10(events) + date + row + col, data = control)
Coefficients:
(Intercept) log10(events) date02/09 date02/11 date02/16
19089.8 3496.2 -6513.6 -4709.4 -6221.7
date02/18 date02/19 date02/20 date02/21 date02/22
-3241.9 -6452.8 -4558.0 -4776.0 -3995.3
date02/23 date03/30 date03/31 rowC rowE
-5054.2 -8507.9 -10061.8 -1882.5 -2429.8
rowG col5 col9
-2641.7 -954.5 -1663.5
rfp.model.1 <- lm(YL2.H ~ log10(events) + date + row + col, data = control)
Model for PHO5pr::RFP
tmp <- dat %>%
# remove one sample with only one valid day of experiment
filter(!(plasmid == "218" & host == "PHO2"), !plasmid %in% c("188", "194", NA)) %>%
nest(data = c(date, BL1.H, YL2.H), .by = c(plasmid, host))
day.var.gfp <- tmp %>%
mutate(model = map(data, function(df) lm(BL1.H ~ date, data = df)),
tidied = map(model, broom::tidy)) %>%
unnest(tidied) %>%
filter(term != "(Intercept)") %>%
mutate(p.adj = p.adjust(p.value, method = "BH")) %>%
select(-data, -model) %>%
filter(p.adj < 0.10) %>%
arrange(plasmid, host)
day.var.rfp <- tmp %>%
mutate(model = map(data, function(df) lm(YL2.H ~ date, data = df)),
tidied = map(model, broom::tidy)) %>%
unnest(tidied) %>%
filter(term != "(Intercept)") %>%
mutate(p.adj = p.adjust(p.value, method = "BH")) %>%
select(-data, -model) %>%
filter(p.adj < 0.10) %>%
arrange(plasmid, host)
there are more systematic shifts in the RFP, significant for row,
col, date and also # of events however, I won’t be removing these
effects yet, because I’ve found that RFP/GFP ratios are pretty
consistent across days. In other words, the variation in GFP and RFP may
be cancelled out.
Check for each plasmid how consistent are the measurements between
days
# extract ximera names
refs <- c("188","194")
# make a test set
day.var.gfp.list <- unique(day.var.gfp$plasmid)
day.var.rfp.list <- unique(day.var.rfp$plasmid)
p <- my_plot_ratio(c(refs,day.var.gfp.list))# +
Warning: There was 1 warning in `mutate()`.
ℹ In argument: `host = fct_recode(host, pho2 = "pho2∆")`.
Caused by warning:
! Unknown levels in `f`: pho2∆
p

High day-to-day GFP variance: 212, 222, 229, 251, 277, 301, 326, 328,
334 High day-to-day RFP variance: 212, 216, 239, 241
Plotting components for chimeras with high day-to-day variance in
Pho4-mNeon
p <- my_plot_ratio(c(refs,day.var.rfp.list))
Warning: There was 1 warning in `mutate()`.
ℹ In argument: `host = fct_recode(host, pho2 = "pho2∆")`.
Caused by warning:
! Unknown levels in `f`: pho2∆
p

Watch out for CSCscC, SCCsS, SCCsS
Plotting components for chimeras with high day-to-day variance in
PHO5pr-mCherry
my_scatter_plot_fix <- function(){
# this function is the same as the one in the script file, but is used to
# plot region 4 effects alone, and doesn't take any input
s1 = my_data_select(pattern = "XXXCCX")
s2 = my_data_select(pattern = "XXXSSX")
scatter.colors = c("ScPho4" = "forestgreen", "CgPho4" = "blue3",
"P2ID:Cg" = "deepskyblue", "P2ID:Sc" = "palegreen2",
"P2ID:mixed" = "gray20")
scatter.size = c("ScPho4" = 3.5, "CgPho4" = 3.5,
"P2ID:Cg" = 2.5, "P2ID:Sc" = 2.5, "P2ID:mixed" = 2.5)
p <- ximera %>%
# exclude the alternative break point sets "A" and "B"
# in particular, pH294 is an alternative break point CCCcs, where P2ID:Cg
# extends to aa 270 instead of 458. it has nearly the same activities as
# CgPho4. However, the CCCCS in the main set has significantly reduced
# activities both with and without Pho4. We later tested whether the
# additional P2ID:Cg4 rescues the effect (see Cg4ext below) and it didn't
# so far, this seems to be an one-off. we need to further investigate its
# activities.
filter(set %in% c("M", "S")) %>%
mutate(A_PHO2 = signif(A_PHO2, digits = 2),
A_pho2 = signif(A_pho2, digits = 2),
group = case_when(
symbol == "CCCCC" ~ "CgPho4",
symbol == "SSSSS" ~ "ScPho4",
plasmid %in% s1 ~ "P2ID:Cg",
plasmid %in% s2 ~ "P2ID:Sc",
.default = "P2ID:mixed"
),
group = fct_relevel(group, names(scatter.colors))) %>%
ggplot(aes(x = A_PHO2, y = A_pho2, label = symbol)) +
geom_abline(slope = 1) +
geom_point(aes(color = group, size = group)) +
scale_color_manual(NULL, values = scatter.colors) +
scale_size_manual(values = scatter.size, guide = "none") +
labs(x = bquote(A[PHO2]), y = bquote(A[pho2*Delta])) +
theme_cowplot() + panel_border(color = "gray30", size = 1.2) +
theme(legend.text = element_text(size = rel(0.8)),
legend.position = c(0.03, 0.83),
axis.title = element_text(face = 2, size = rel(1.2)),
axis.line = element_blank())
return(p)
}
most of the day-to-day variance are canceled out after RFP/GFP
normalization
All chimera, scatter plot
Plot all chimeras, coloring based on P2ID source
p <- my_scatter_plot_fix()
ggsave(filename = "../img/20240308-all-chimera-scatter-color-by-P2ID.png",
plot = p, width = 4.5, height = 4, dpi = 300)
ggplotly(p + labs(x = "A<sub>PHO2</sub>", y = "A<sub>pho2</sub>") +
theme_gray(base_size = 16) +
theme(legend.text = element_markdown()),
tooltip = c("label", "x", "y"))
my_scatter_plot_all <- function(){
# this function is the same as my_scatter_plot_fix except that it plots all the chimeras
# without coloring them differently. for figure 5
s1 = my_data_select(pattern = "XXXCCX")
s2 = my_data_select(pattern = "XXXSSX")
scatter.colors = c("ScPho4" = "forestgreen", "CgPho4" = "blue3",
"P2ID:Cg" = "gray20", "P2ID:Sc" = "gray20",
"P2ID:mixed" = "gray20")
scatter.size = c("ScPho4" = 3.5, "CgPho4" = 3.5,
"P2ID:Cg" = 2.5, "P2ID:Sc" = 2.5, "P2ID:mixed" = 2.5)
p <- ximera %>%
filter(set %in% c("M", "S")) %>%
mutate(A_PHO2 = signif(A_PHO2, digits = 2),
A_pho2 = signif(A_pho2, digits = 2),
group = case_when(
symbol == "CCCCC" ~ "CgPho4",
symbol == "SSSSS" ~ "ScPho4",
plasmid %in% s1 ~ "P2ID:Cg",
plasmid %in% s2 ~ "P2ID:Sc",
.default = "P2ID:mixed"
),
group = fct_relevel(group, names(scatter.colors))) %>%
ggplot(aes(x = A_PHO2, y = A_pho2, label = symbol)) +
geom_abline(slope = 1) +
geom_point(aes(color = group, size = group)) +
scale_color_manual(NULL, values = scatter.colors) +
scale_size_manual(values = scatter.size, guide = "none") +
labs(x = bquote(A[PHO2]), y = bquote(A[pho2])) +
theme_cowplot() + panel_border(color = "gray30", size = 1.2) +
theme(legend.text = element_text(size = rel(0.8)),
legend.position = "none",
axis.title = element_text(face = 2, size = rel(1.2)),
axis.line = element_blank())
return(p)
}
this function is the same as my_scatter_plot_fix except that it plots
all the chimeras without coloring them differently. for figure 5
p <- my_scatter_plot_all()
ggsave(filename = "../img/20241121-all-chimera-scatter.png",
plot = p, width = 4.5, height = 4, dpi = 300)
ggplotly(p + labs(x = "A<sub>PHO2</sub>", y = "A<sub>pho2</sub>") +
theme_gray(base_size = 16) +
theme(legend.text = element_markdown()),
tooltip = c("label", "x", "y"))
Plot all chimeras, for Fig. 5
my_plot_subset_ximera <- function(symbols){
# this function plots a subset of the chimeras as horizontal bar plots
# showing the Rel. A_PHO2 and %A_pho2∆ values
# it takes as input a vector containing the symbols for the chimeras for
# plotting. the order in the vector determines the plot order
# the endogenous ScPho4 and CgPho4 are implied
missing <- setdiff(symbols, ximera$symbol)
if(length(missing) != 0)
stop(paste(missing, "are not found", sep = " "))
tmp <- filter(ximera, symbol %in% c("SSSSS", "CCCCC", symbols)) %>%
mutate(
rSE_PHO2 = se_PHO2 / A_PHO2[symbol == "SSSSS"],
rSE_pho2 = se_pho2 / A_pho2[symbol == "SSSSS"]
) %>%
pivot_longer(cols = c(rA_PHO2, rA_pho2, rSE_PHO2, rSE_pho2),
#pivot_longer(cols = c(A_PHO2, A_pho2, se_PHO2, se_pho2),
names_to = c(".value", "parameter"), names_sep = "_",
values_to = "value") %>%
mutate(parameter = fct_relevel(parameter, "PHO2"),
symbol = factor(symbol, levels =
unique(c("SSSSS", "CCCCC", symbols)))) %>%
select(-c(A_PHO2:boost))
# labeller
par.explain <- c(
PHO2 = "Rel. A<sub>PHO2</sub>",
#boost = "Boost",
pho2 = "Rel. A<sub>pho2∆</sub>"
)
p <- ggplot(tmp, aes(y = symbol, x = rA)) +
geom_col(width = 0.5, color = "black", fill = "gray80") +
geom_vline(xintercept = 1, linetype = 2, color = "gray30") +
geom_errorbar(aes(xmin = rA - rSE, xmax = rA + rSE), width = 0.2) +
facet_wrap(~parameter, scales = "free_x",# switch = "x",
labeller = labeller(parameter = par.explain)) +
scale_y_discrete(limits = rev) +
scale_x_continuous(expand = expansion(mult = c(0.02, 0.05))) +
theme_cowplot() + panel_border(color = "gray30") +
background_grid(major = "y", minor = "none") +
theme(axis.text.y = element_text(family = "courier"),
axis.title = element_blank(),
axis.line = element_blank(),
strip.placement = "outside",
strip.background = element_blank(),
strip.text = element_markdown())
return(p)
}
Spotlight individual chimeras
The goal here is to plot individual chimeras in order to test
specific hypotheses and make certain points.
- We separately tested and found that CgPho4 DBD binds the consensus
DNA more strongly than ScPho4 does, and it also has two additional
activation booster regions, which enhance the activity of the main AD.
We therefore hypothesize that by replacing the corresponding regions in
ScPho4 with the parts from CgPho4, we would create a chimeric TF that is
not or far less dependent on Pho2.
- We also expect that those regions additively contribute to the
reduced Pho2-dependence, shown as increased TF activity of the chimera
in the pho2∆ background.
Design plot
my_plot_subset_ximera_alt <- function(symbols){
# this function plots a subset of the chimeras as horizontal bar plots
# showing the Rel. A_PHO2 and %A_pho2∆ values
# it takes as input a vector containing the symbols for the chimeras for
# plotting. the order in the vector determines the plot order
# the endogenous ScPho4 and CgPho4 are implied
missing <- setdiff(symbols, ximera$symbol)
if(length(missing) != 0)
stop(paste(missing, "are not found", sep = " "))
tmp <- filter(dat_sep, symbol %in% c("SSSSS", "CCCCC", symbols)) %>%
mutate(host = fct_relevel(host, "PHO2"),
symbol = factor(symbol, levels =
unique(c("SSSSS", "CCCCC", symbols))))
tmp %>% count(symbol, host) %>% print()
# labeller
par.explain <- c(
PHO2 = "A<sub>PHO2</sub>",
#boost = "Boost",
pho2 = "A<sub>pho2∆</sub>"
)
p <- ggplot(tmp, aes(y = symbol, x = A)) +
geom_bar(stat = "summary", fun = "mean",
width = 0.5, color = "black", fill = "gray80") +
stat_summary(fun.data = "mean_cl_boot", geom = "linerange",
color = "steelblue4") +
geom_point(data = filter(tmp, !symbol %in% c("CCCCC", "SSSSS")),
size = 0.6, shape = 3, color = "gray30") +
#geom_vline(xintercept = 1, linetype = 2, color = "gray30") +
#geom_errorbar(aes(xmin = rA - rSE, xmax = rA + rSE), width = 0.2) +
facet_wrap(~host, scales = "free_x",# switch = "x",
labeller = labeller(host = par.explain)) +
scale_y_discrete(limits = rev) +
scale_x_continuous(expand = expansion(mult = c(0.02, 0.05))) +
theme_cowplot() + panel_border(color = "gray30") +
background_grid(major = "y", minor = "none") +
theme(axis.text.y = element_text(family = "courier"),
axis.title = element_blank(),
axis.line = element_blank(),
strip.placement = "outside",
strip.background = element_blank(),
strip.text = element_markdown())
return(p)
}
Update 2024-11-22
Alternative design, with individual points and not relative to
ScPho4. Also, individual datapoints were plotted for samples with <10
replicates
selected <- as.character(
expression(CSSSS, SCSSS, CCSSS, SSCSS, CSCSS, SCCSS, CCCSS, SSSSC, SSCSC, CSCSC, CSScC))
#selected <- filter(meta, symbol %in% selected) %>% pull(plasmid)
my_plot_subset_ximera(selected)
ggsave("../img/20240308-selected-chimera-rel-activity.png", width = 4, height = 4)

# plot with absolute A not relative, and plot individual data points
my_plot_subset_ximera_alt(selected)
ggsave("../img/20241122-selected-chimera-rel-activity.png", width = 4, height = 4)

Minimal CgPho4 parts for A_pho2
The chimera with the least amount of CgPho4 and yet have appreciable
activity in the absence of Pho2 is These include SSSSS, CCCCC, SSSSC,
CSSSS, SSCSS, CSCSS, CSSSC, CSCSC
# select the chimeras
selected <- as.character(
expression(SSSSS, CSSSS, SCSSS, SSCSS, CCSSS, CSCSS, SCCSS, CCCSS)
)
# extract the data
tmp <- ximera %>%
filter(symbol %in% selected, set == "M") %>%
select(plasmid, symbol, group) %>%
inner_join(dat, by = "plasmid") %>%
mutate( `R/G` = YL2.H / BL1.H ) %>%
filter(flag == "pass") %>%
select(-nRFP, -nGFP, -well, -flag)
# prepare the factor levels
split <- c(1,1,1,2); names(split) <- c("R1", "AD", "NLS")
tmp <- tmp %>%
separate_wider_position(symbol, split) %>%
mutate(across(R1:NLS, ~factor(.x, levels = c("S", "C"))))
# test A_PHO2
print("Testing A_PHO2")
[1] "Testing A_PHO2"
lm.res <- tmp %>%
filter(host == "PHO2") %>%
lm(`R/G` ~ (R1*AD*NLS), data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ (R1 * AD * NLS), data = .)
Residuals:
Min 1Q Median 3Q Max
-6.4865 -1.6014 -0.5069 1.0472 12.7647
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 1.764e+01 2.391e-01 7.375e+01 5.546e-137 0.000
R1C 3.591e+00 1.196e+00 3.004e+00 3.046e-03 0.012
ADC 2.237e+00 1.196e+00 1.871e+00 6.296e-02 0.063
NLSC 1.499e+01 1.196e+00 1.254e+01 2.350e-26 0.000
R1C:ADC 7.578e+00 2.109e+00 3.593e+00 4.213e-04 0.002
R1C:NLSC -5.110e+00 1.904e+00 -2.684e+00 7.959e-03 0.021
ADC:NLSC 5.586e+00 2.043e+00 2.734e+00 6.873e-03 0.021
R1C:ADC:NLSC -1.383e+01 3.064e+00 -4.514e+00 1.142e-05 0.000
Residual standard error: 2.869 on 181 degrees of freedom
Multiple R-squared: 0.825, Adjusted R-squared: 0.8182
F-statistic: 121.9 on 7 and 181 DF, p-value: < 2.2e-16
# store the test results for plotting
res.PHO2 <- coef(lm.res)[-1,] %>% as_tibble(rownames = "component")
# test A_pho2
print("Testing A_pho2∆")
[1] "Testing A_pho2∆"
lm.res <- tmp %>%
filter(host == "pho2") %>%
lm(`R/G` ~ (R1*AD*NLS), data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ (R1 * AD * NLS), data = .)
Residuals:
Min 1Q Median 3Q Max
-0.63622 -0.21548 -0.01744 0.10489 1.15930
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 1.072e+00 5.581e-02 1.920e+01 1.312e-29 0.00
R1C 6.167e-01 1.477e-01 4.176e+00 8.404e-05 0.00
ADC -1.153e-02 1.477e-01 -7.805e-02 9.380e-01 1.00
NLSC 1.819e+00 1.477e-01 1.232e+01 3.389e-19 0.00
R1C:ADC 5.923e-01 2.433e-01 2.435e+00 1.746e-02 0.07
R1C:NLSC 1.479e+00 2.433e-01 6.081e+00 5.623e-08 0.00
ADC:NLSC 1.575e-01 2.433e-01 6.475e-01 5.194e-01 1.00
R1C:ADC:NLSC -3.795e-01 3.660e-01 -1.037e+00 3.033e-01 0.91
Residual standard error: 0.3349 on 70 degrees of freedom
Multiple R-squared: 0.9553, Adjusted R-squared: 0.9509
F-statistic: 213.9 on 7 and 70 DF, p-value: < 2.2e-16
# store the test results for plotting
res.pho2 <- coef(lm.res)[-1,] %>% as_tibble(rownames = "component")
# combine the results
test.res <- bind_rows(
"A_PHO2" = res.PHO2, "A_pho2" = res.pho2, .id = "parameter"
)
Region 1-3 main effects and interactions
par.explain <- c(
A_PHO2 = "A<sub>PHO2</sub>",
A_pho2 = "A<sub>pho2∆</sub>"
)
p <- test.res %>%
rename(estimate = Estimate, se = `Std. Error`) %>%
mutate(
component = gsub("C", "", component) %>% fct_inorder(),
parameter = factor(parameter, levels = c("A_PHO2", "A_pho2")),
sig = P.adj < 0.05
) %>%
ggplot(aes(x = component, y = estimate)) +
geom_hline(yintercept = 0, linetype = 1, color = "gray50") +
geom_col(aes(fill = P.adj < 0.05), width = 0.5, color = "black") +
geom_pointrange(aes(ymin = estimate-se, ymax = estimate+se), size = 0.2) +
facet_wrap(~parameter, scales = "free_y", nrow = 2,
labeller = labeller(parameter = par.explain)) +
scale_x_discrete() +
scale_y_continuous() +
scale_fill_manual(NULL,
values = c("gray90", "gray50")) +
theme_cowplot() + panel_border(color = "gray30") +
background_grid(major = "y", minor = "y") +
theme(axis.text.x = element_text(angle = 90, hjust = 1, size = rel(1)),
axis.title = element_blank(),
axis.line = element_blank(),
legend.position = "bottom",
strip.placement = "outside",
strip.background = element_blank(),
strip.text = element_markdown(size = rel(1)))
p
ggsave("../img/20240404-region1-3-epistasis-plot.png", width = 3.5, height = 4.5)

Plot the result
selected <- as.character(
expression(CCCSC, CCCcsC, CCCscC,
CSSCC, CSSSC, CSScsC, CSSscC
))
#selected <- filter(meta, symbol %in% selected) %>% pull(plasmid)
my_plot_subset_ximera_alt(selected) + scale_x_continuous(expand = expansion(mult = c(0.02, 0.15)))
Scale for x is already present.
Adding another scale for x, which will replace the existing scale.
ggsave("../img/20241122-P2ID-split-rel-activity.png", width = 3.5, height = 3.5)

Region 4 splits
We have so far focused on the main set with the 5 region design. In
the scatter plot below, we see that there is a subset of chimeras in
between the P2ID:Sc and P2ID:Cg ones. They are interesting in that their
A_pho2∆/A_PHO2 ratios are intermediate.
scatter
# select the chimeras
selected <- as.character(
expression(CCCCC, CCCSC, CCCcsC, CCCscC)
)
# extract the data
tmp <- ximera %>%
filter(symbol %in% selected) %>%
select(plasmid, symbol, group) %>%
inner_join(dat, by = "plasmid") %>%
mutate( `R/G` = YL2.H / BL1.H ) %>%
filter(flag == "pass") %>%
select(-nRFP, -nGFP, -well, -flag) %>%
mutate(symbol = factor(symbol, levels = !!selected))
# test A_PHO2
print("Testing A_PHO2")
[1] "Testing A_PHO2"
lm.res <- tmp %>%
filter(host == "PHO2") %>%
lm(`R/G` ~ symbol, data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ symbol, data = .)
Residuals:
Min 1Q Median 3Q Max
-5.2250 -1.0489 -0.0647 1.2813 6.6290
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 1.809e+01 3.445e-01 5.252e+01 2.067e-45 0
symbolCCCSC 4.207e+00 9.115e-01 4.615e+00 2.766e-05 0
symbolCCCcsC -7.086e+00 9.115e-01 -7.774e+00 3.722e-10 0
symbolCCCscC 1.362e+01 9.115e-01 1.494e+01 4.511e-20 0
Residual standard error: 2.067 on 50 degrees of freedom
Multiple R-squared: 0.8711, Adjusted R-squared: 0.8633
F-statistic: 112.6 on 3 and 50 DF, p-value: < 2.2e-16
# store the test results for plotting
#res.PHO2 <- coef(lm.res)[-1,] %>% as_tibble(rownames = "component")
# test A_pho2
print("Testing A_pho2∆")
[1] "Testing A_pho2∆"
lm.res <- tmp %>%
filter(host == "pho2") %>%
lm(`R/G` ~ symbol, data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ symbol, data = .)
Residuals:
Min 1Q Median 3Q Max
-5.1275 -0.5912 -0.0824 0.5645 5.2165
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 1.674e+01 2.869e-01 5.834e+01 1.176e-47 0.000
symbolCCCSC -1.195e+01 7.591e-01 -1.574e+01 5.226e-21 0.000
symbolCCCcsC -1.255e+01 7.591e-01 -1.654e+01 6.456e-22 0.000
symbolCCCscC 1.928e+00 7.591e-01 2.540e+00 1.424e-02 0.014
Residual standard error: 1.721 on 50 degrees of freedom
Multiple R-squared: 0.9093, Adjusted R-squared: 0.9038
F-statistic: 167 on 3 and 50 DF, p-value: < 2.2e-16
Statistical tests for group 1
# select the chimeras
selected <- as.character(
expression( CSSCC, CSSSC, CSScsC, CSSscC )
)
# extract the data
tmp <- ximera %>%
filter(symbol %in% selected) %>%
select(plasmid, symbol, group) %>%
inner_join(dat, by = "plasmid") %>%
mutate( `R/G` = YL2.H / BL1.H ) %>%
filter(flag == "pass") %>%
select(-nRFP, -nGFP, -well, -flag) %>%
mutate(symbol = factor(symbol, levels = !!selected))
# test A_PHO2
print("Testing A_PHO2")
[1] "Testing A_PHO2"
lm.res <- tmp %>%
filter(host == "PHO2") %>%
lm(`R/G` ~ symbol, data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ symbol, data = .)
Residuals:
Min 1Q Median 3Q Max
-3.4560 -0.5974 -0.1128 0.4629 4.4143
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 6.270e+00 7.393e-01 8.481e+00 4.673e-08 0.000
symbolCSSSC 8.962e+00 1.046e+00 8.572e+00 3.949e-08 0.000
symbolCSScsC 2.558e+00 1.046e+00 2.447e+00 2.378e-02 0.024
symbolCSSscC 6.372e+00 1.046e+00 6.094e+00 5.891e-06 0.000
Residual standard error: 1.811 on 20 degrees of freedom
Multiple R-squared: 0.8127, Adjusted R-squared: 0.7846
F-statistic: 28.93 on 3 and 20 DF, p-value: 1.791e-07
# store the test results for plotting
#res.PHO2 <- coef(lm.res)[-1,] %>% as_tibble(rownames = "component")
# test A_pho2
print("Testing A_pho2∆")
[1] "Testing A_pho2∆"
lm.res <- tmp %>%
filter(host == "pho2") %>%
lm(`R/G` ~ symbol, data = .) %>%
summary()
# adding adjusted P-value
lm.res$coefficients <- cbind(
coef(lm.res),
"P.adj" = p.adjust(coef(lm.res)[,'Pr(>|t|)'], method = "holm")
)
print(lm.res)
Call:
lm(formula = `R/G` ~ symbol, data = .)
Residuals:
Min 1Q Median 3Q Max
-0.77040 -0.08856 -0.03391 0.16099 0.37316
Coefficients:
Estimate Std. Error t value Pr(>|t|) P.adj
(Intercept) 5.244e+00 1.082e-01 4.846e+01 3.257e-22 0.000
symbolCSSSC -2.822e+00 1.530e-01 -1.844e+01 5.058e-14 0.000
symbolCSScsC -3.061e+00 1.530e-01 -2.000e+01 1.076e-14 0.000
symbolCSSscC -1.243e-01 1.530e-01 -8.125e-01 4.261e-01 0.426
Residual standard error: 0.265 on 20 degrees of freedom
Multiple R-squared: 0.9726, Adjusted R-squared: 0.9685
F-statistic: 237.1 on 3 and 20 DF, p-value: 8.565e-16
Statistical tests for group 2
split <- c(1,1,1,1,1); names(split) <- paste0("P", 1:5)
tmp <- ximera %>%
filter(set == "M", group != "n.f.") %>%
separate_wider_position(symbol, split) %>%
mutate(across(P1:P5, ~factor(.x, levels = c("S", "C"))))
lm <- lm(A_pho2 ~ (P1+P2+P3+P4+P5), data = tmp)
summary(lm)
Call:
lm(formula = A_pho2 ~ (P1 + P2 + P3 + P4 + P5), data = tmp)
Residuals:
Min 1Q Median 3Q Max
-3.5062 -1.1430 -0.0082 0.8638 6.5285
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.7264 1.0811 -0.672 0.5102
P1C 2.3316 0.8687 2.684 0.0152 *
P2C 0.7880 0.8687 0.907 0.3763
P3C 2.5260 0.9010 2.804 0.0117 *
P4C 4.0647 0.9792 4.151 0.0006 ***
P5C 1.2254 0.9010 1.360 0.1906
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2.086 on 18 degrees of freedom
Multiple R-squared: 0.7398, Adjusted R-squared: 0.6675
F-statistic: 10.23 on 5 and 18 DF, p-value: 9.054e-05
Region main effect
split <- c(1,1,1,1,1); names(split) <- paste0("P", 1:5)
tmp <- ximera %>%
filter(set == "M", group != "n.f.") %>%
separate_wider_position(symbol, split) %>%
mutate(across(P1:P5, ~factor(.x, levels = c("S", "C"))))
lm <- lm(A_pho2 ~ (P1+P2+P3+P4+P5), data = tmp)
summary(lm)
The main effects were calculated by averaging over all chimeras with
CgPho4 region at the respective position. I’d like to break them down by
backgrounds. For example, for region 3, I’d like to see the pairwise
comparisons between CCCSS and CCSSS, where only region 3 differs. The
steps are
- select the region to be compared. split the symbol into two parts -
the genotype of the focal region and the rest
- group by the second part (rest) and calculate the differential
my_calc_region_effect <- function(region, variable){
# this function takes the name of a variable of interest
# x specifies the foreground region, which will be examined for its effect on
# the variable of interest.
# it then transforms the ximera data frame to preserve only the variable of
# interest, pivots it wider after grouping by the background composition.
# prepare the data by mutating the symbol column into fg and bg
valid.var <- c("A_PHO2", "A_pho2", "rA_PHO2", "rA_pho2", "boost")
if(!variable %in% valid.var)
stop(paste0("Please specify one of the valid variable names:",
paste(valid.var, collapse = ", ")))
tmp <- ximera %>%
filter(set == "M") %>%
select(plasmid, symbol, var = {{ variable }}) %>%
mutate(fg = str_sub(symbol, region, region) %>% toupper(),
bg = symbol %>% toupper())
# replace the foreground region with X for grouping
str_sub(tmp$bg, region, region) <- "X"
# reorganize the tibble for easier handling, optional
tmp <- relocate(tmp, fg, bg, .before = symbol) %>% select(-symbol)
# pivot the data into a wide format such that for each background, there
# are two values for the variable of interest, one from the chimera with
# CgPho4's version in the foreground and another with ScPho4's version
tmp <- tmp %>%
select(plasmid, fg, bg, var) %>%
pivot_wider(id_cols = bg, names_from = "fg",
values_from = c(plasmid, var)) %>%
unite(plasmid, starts_with("plasmid")) %>%
mutate(label = paste(bg, plasmid, sep = "\n"))
return(tmp)
}
x = 5
p1 <- my_plot_region_effect_onevar(x, "A_PHO2")
p2 <- my_plot_region_effect_onevar(x, "A_pho2")
subplot(p1, p2, margin = 0.05) %>%
layout(title = paste("Region", x, "swap effect on A_PHO2 and A_pho2", sep = " "),
xaxis = list(title = paste0("Region ", x, " from CgPho4")),
yaxis = list(title = paste0("Region ", x, " from ScPho4")) )
Here, I’d like to take what I build above and create a new tibble, in
which each row is a different background (makeup of the chimera except
for the focal region). The value columns are:
- dA_PHO2 = A_PHO2_Cg - A_PHO2_Sc
- dA_pho2 = A_pho2_Cg - A_pho2_Sc
- A_PHO2_Sc = A_PHO2_Sc
The goal is to plot dA_PHO2 and dA_pho2 side-by-side for each
background.
my_comp_region_effect <- function(region){
# this function uses my_calc_region_effect to get the value for the variable of interest
# with either Cg or Sc version in the focal region, separately for each background composition
# it does so for two variables, A_PHO2 and A_pho2, then calculate dA_PHO2, dA_pho2, and
# combine them
PHO2 = my_calc_region_effect(region, "A_PHO2") %>%
mutate(dA_PHO2 = var_C - var_S,
# mean A_PHO2
M_PHO2 = (var_S + var_C)/2,
NF = ifelse(M_PHO2 <=3.5, TRUE, FALSE)) %>%
select(-var_S, -var_C)
pho2 = my_calc_region_effect(region, "A_pho2") %>%
mutate(dA_pho2 = var_C - var_S,
M_pho2 = (var_S + var_C)/2) %>%
select(-var_S, -var_C)
dat <- full_join(PHO2, pho2, by = c("bg", "plasmid", "label")) %>%
select(bg, plasmid, dA_PHO2, dA_pho2, M_PHO2, M_pho2, NF)
return(dat)
}
my_plot_region_effect_twovar_line("1", "4")# %>% ggplotly()
ggsave("../img/20240310-region-swap-effect-1-on-4.png", width = 6, height = 4)
my_plot_region_effect_twovar_line("3", "4")# %>% ggplotly()
ggsave("../img/20240310-region-swap-effect-3-on-4.png", width = 6, height = 4)
my_plot_region_effect_twovar_line("4", "5")# %>% ggplotly()
ggsave("../img/20231221-region-swap-effect-4-on-5.png", width = 6, height = 4, dpi = 150)
my_plot_region_effect_twovar_line("5", "4")# %>% ggplotly()
ggsave("../img/20231224-region-swap-effect-5-on-4.png", width = 6, height = 4, dpi = 200)
The main plotting functions are now in a separate script file in
../script. The plotting function below is to adapt the plot
for a figure in the paper, simultaneously showing regions 1-3.
my_plot_region_effect_twovar_line_par <- function(regions){
# this function uses my_comp_region_effect to generate the data
# and plot the difference in A_PHO2 and A_pho2 between the CgPho4 vs ScPho4
# in the focal region
dat <- map_dfr(regions, \(region) my_comp_region_effect(region), .id = "region") %>%
pivot_longer(cols = c(dA_PHO2, dA_pho2),
names_to = "host", values_to = "diff") %>%
mutate(host = fct_recode(host, `PHO2` = "dA_PHO2", `pho2∆` = "dA_pho2"),
host = fct_relevel(host, "PHO2"))
# specify grouping variable
dat <- mutate(dat,
grp = str_sub(bg, 4, 4) %>% toupper(),
grp = fct_recode(grp, CgPho4 = "C", ScPho4 = "S"))#,
#sh = str_sub(bg, 5, 5) %>% toupper(),
#sh = fct_recode(sh, CgPho4 = "C", ScPho4 = "S") )
# specify arrow annotation
arrow.x = 0.7
arrow.y = (max(dat$diff) - min(dat$diff)) / 5
# plot
p <- dat %>%
ggplot(aes(x = host, y = diff, label = bg)) +
geom_point(aes(color = grp), size = 2, alpha = 0.8,
position = position_jitter(0.1)) +
geom_line(aes(group = bg), linewidth = 0.2, alpha = 0.8) +
facet_grid(region ~ grp, labeller = labeller(
grp = c(CgPho4 = "P2ID:Cg", ScPho4 = "P2ID:Sc"),
region = label_both
)) +
scale_color_manual("P2ID:", values = c("orange", "gray30"), guide = "none") +
#scale_shape_manual("DBD:", values = c(19, 1)) +
ylab("Region swap effect (Cg-Sc)") +
theme_bw(base_size = 18) +
theme(
axis.title.x = element_blank(),
axis.title.y = element_text(size = rel(0.9)),
axis.text.x = element_text(face = 3),
axis.text.y = element_text(size = rel(0.8)),
legend.text = element_text(size = rel(0.8)),
legend.title = element_text(size = rel(0.9)),
legend.position = "top",
strip.background = element_blank()
)
return(p)
}
my_plot_region_effect_twovar_line_par(c(1,3))
ggsave("../img/20240310-region-swap-effect-1n3-on-4.png",
width = 5, height = 3.5)
my_scatter_plot("XXXXCS") +
labs(x = bquote(A[PHO2]), y = bquote(A[pho2]))
ggsave(filename = "../img/20240214-Pho4-chimeras-scatter-P2ID_Cg-DBD_Sc.png",
plot = p, width = 6, height = 4, dpi = 300)

P2ID:Cg_DBD:Sc fail
Highlight the subset of the chimeras with P2ID:Cg + DBD:Sc, most of
which are non functional
my_scatter_plot("XXXXCS") +
labs(x = bquote(A[PHO2]), y = bquote(A[pho2]))
ggsave(filename = "../img/20240214-Pho4-chimeras-scatter-P2ID_Cg-DBD_Sc.png",
plot = p, width = 6, height = 4, dpi = 300)
x <- my_data_select(pattern = "XXXXCS", Set = "M")
my_data_prep(x) %>%
mutate(group = fct_recode(group, "chimera" = "n.f.")) %>%
my_plot_components()
ggsave("../img/20240213-P2ID_Cg-DBD_Sc-components.png", width = 8, height = 5)
Triangle heatmap
First, write a function to generate the data for plotting. If we are
going to use ggplot, we need a tibble to store the data, something in
the following form
| 209 |
CCSCC |
3 |
3 |
8.25 |
7.82 |
0.468 |
1.06 |
0.94 |
If we are ok with using non ggplot - heatmaps are not ggplot’s
strength anyways - we can just build a matrix.
Note that this way of summarizing the data has many limitaitons: 1)
it requires specifying the reference, either CCCCC or SSSSS. Everything
is measured against that; 2) it only shows pairwise (two region)
interactions. This turns out to be fine with five regions, since every
chimera can be expressed as either a 0, 1 or 2 region swap from one of
the two reference genotypes. With 6 or more regions, higher level (3 or
more region) interactions cannot be visualized this way. Because of
this, we will focus on just the main set for this analysis.
To build the matrix, we need to first identify the chimeras that
belong to the set. For that, we will use the “main” set, with the five
region split, for the moment at least. The function will first determine
which reference to use. If we use SSSSS as the reference, for example,
we will assign 0 to the reference. All other chimeras with 1 or 2
regions from Cg will be used to fill an upper triangular matrix, using
one of the values of interest, e.g., A_PHO2.
my_upper_triangular_mat <- function(alt = "C", var = "A_PHO2", nf.as.na = F){
# given the alternative allele (C/S) and a variable of interest, e.g., A_PHO2,
# output an upper triangular matrix containing the values from the variable
# of interest, with the row and col numbers based on the first and second
# positions containing the alternative allele. If all positions contain the
# reference allele, the value is subtracted from all values in the matrix
# when just one position is the alternative allele, the value in the diagonal
# is set. when there are more than 2 regions containing the alternative allele
# skip.
# if "nf.as.na = TRUE", evaluate if the activity of either of the two chimeras
# being compared is non functional. if yes, set the corresponding matrix value
# to NA
out_mat <- matrix(NA, nrow = 5, ncol = 5)
ref_val <- NA
dat <- filter(ximera, set == "M") %>%
mutate(S = as.character(symbol) %>% toupper())
if(nf.as.na){
dat <- filter(dat, group != "n.f.")
}
for(i in seq(1, nrow(dat))){
symbol = dat[i, "S"]
# determine which positions contain the alternative allele
p = str_locate_all(symbol, alt)[[1]][,"start"]
l = length(p) # how many positions contain the alt allele
v = dat[[var]][i] # retrieve the value of the variable
if(l == 0)
ref_val = v
else if(l == 1)
out_mat[p, p] = v
else if(l == 2)
out_mat[p[1], p[2]] = v
}
out_mat = out_mat - ref_val
return(out_mat)
}
my_combined_triangular_mat <- function(alt = "C"){
# given the alternative allele (C/S), output a matrix containing the values
# for both with and without Pho2, arranged in two complementary triagular
# matrices, with the row and col numbers based on the first and second
# positions containing the alternative allele. If all positions contain the
# reference allele, the value is subtracted from all values in the matrix
# when just one position is the alternative allele, the value in the diagonal
# is set. when there are more than 2 regions containing the alternative allele
# skip.
out_mat <- matrix(NA, nrow = 6, ncol = 6)
upper <- cbind(NA, my_upper_triangular_mat(alt, var = "A_PHO2", )) %>%
rbind(., NA)
lower <- rbind(NA, t(my_upper_triangular_mat(alt, var = "A_pho2"))) %>%
cbind(., NA)
out_mat = ifelse(is.na(upper), lower, upper)
return(out_mat)
}
my_plot_triangle_heatmap <- function(alt, var){
# this function takes the output of the function above and makes a heatmap
# using pheatmap function, then rotates it using grid graphics
# thanks to https://bookdown.org/rdpeng/RProgDA/the-grid-package.html#grid-graphics-coordinate-systems
# adding title based on https://davetang.github.io/muse/pheatmap.html
# construct title of plot
ref = ifelse(alt == "C", "ScPho4", "CgPho4")
bg = ifelse(var == "A_PHO2", "with PHO2", "w/o pho2")
my_title <- paste("Epistasis between regions on", ref, "background", bg)
test <- my_upper_triangular_mat(alt = alt, var = var)
paletteLength = 50
myColors <- colorRampPalette(c("steelblue3", "gray90", "red"))(paletteLength)
rng <- max(abs(test), na.rm = TRUE)
myBreaks <- c(seq(-rng, 0, length.out=ceiling(paletteLength/2) + 1),
seq(rng/paletteLength, rng,
length.out=floor(paletteLength/2)))
p <- pheatmap::pheatmap(test, color = myColors, breaks = myBreaks,
border_color = NA, na_col = NA, silent = TRUE,
cluster_cols = FALSE, cluster_rows = FALSE)
vp <- viewport(x = 0.5, y = 0.25,
width = unit(4.5, "in"), height = unit(4.5, "in"), angle = 47)
grid.newpage()
pushViewport(vp)
grid.draw(p$gtable)
popViewport()
grid.text(label = my_title, x = 0.5, y = 0.95, gp = gpar(fontsize = 16, fontface = "bold"))
return(p)
}
my_plot_combined_triangle_heatmap <- function(alt){
# this function takes the output of the function my_combined_triangular_mat()
# using pheatmap function, then rotates it using grid graphics
# thanks to https://bookdown.org/rdpeng/RProgDA/the-grid-package.html#grid-graphics-coordinate-systems
# adding title based on https://davetang.github.io/muse/pheatmap.html
# construct title of plot
ref = ifelse(alt == "C", "ScPho4", "CgPho4")
my_title <- paste("Epistasis between regions on", ref, "background")
test <- my_combined_triangular_mat(alt = alt)
paletteLength = 50
myColors <- colorRampPalette(c("steelblue", "gray90", "red"))(paletteLength)
rng <- max(abs(test), na.rm = TRUE)
myBreaks <- c(seq(-rng, 0, length.out=ceiling(paletteLength/2) + 1),
seq(rng/paletteLength, rng,
length.out=floor(paletteLength/2)))
p <- pheatmap::pheatmap(test, color = myColors, breaks = myBreaks,
border_color = NA, na_col = NA, silent = TRUE,
cluster_cols = FALSE, cluster_rows = FALSE)
vp <- viewport(x = 0.5, y = 0.45,
width = unit(3, "in"), height = unit(2.8, "in"), angle = 47)
grid.newpage()
pushViewport(vp)
grid.draw(p$gtable)
popViewport()
grid.text(label = my_title, x = 0.5, y = 0.95,
gp = gpar(fontsize = 16, fontface = "bold"))
grid.text(label = "With Pho2", x = 0.1, y = 0.65, just = c("left", "top"),
gp = gpar(fontsize = 14, fontface = "bold"))
grid.text(label = "Without pho2", x = 0.1, y = 0.25, just = c("left", "top"),
gp = gpar(fontsize = 14, fontface = "bold"))
return(p)
}
png("../img/20240115-triangle-heatmap-CgPho4-ref.png", width = 7, height = 5, units = "in", res = 300)
p1 <- my_plot_combined_triangle_heatmap("C")
dev.off()
png("../img/20240115-triangle-heatmap-ScPho4-ref.png", width = 7, height = 5, units = "in", res = 300)
p2 <- my_plot_combined_triangle_heatmap("S")
dev.off()
LS0tCnRpdGxlOiAiRTAxMyBQaG80IGNoaW1lcmEgYWN0aXZpdHkgYW5hbHlzaXMgdXNpbmcgUEhPNSByZXBvcnRlciwgYW5hbHlzaXMiCmF1dGhvcjogIkJpbiBIZSIKZGF0ZTogIjIwMjMtMTAtMzEgdXBkYXRlZCBgciBTeXMuRGF0ZSgpYCIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHllcwogICAgdG9jX2Zsb2F0OiB5ZXMKICAgIGNvZGVfZm9sZGluZzogaGlkZQogIHBkZl9kb2N1bWVudDoKICAgIHRvYzogeWVzCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogeWVzCiAgICBkZl9wcmludDogcGFnZWQKLS0tCgpgYGB7ciBtZXNzYWdlPUZBTFNFfQpyZXF1aXJlKHBsb3RseSkKcmVxdWlyZSh0aWR5dmVyc2UpCnJlcXVpcmUoZ2dyaWRnZXMpCnJlcXVpcmUoY293cGxvdCkKcmVxdWlyZShSQ29sb3JCcmV3ZXIpCnJlcXVpcmUoZ3JpZCkKcmVxdWlyZShnZ3RleHQpCmBgYAoKYGBge3J9Cm9sZCA8LSB0aGVtZV9zZXQodGhlbWVfYncoYmFzZV9zaXplID0gMTYpKQpgYGAKClRoaXMgaXMgdGhlIHNlY29uZCBwYXJ0IG9mIHRoZSBhbmFseXNpcy4gSW4gdGhlIGZpcnN0IHBhcnQgKHNlZSBgLi4vaW5wdXQvUEhPNS1kYXRhLzIwMjMxMDE5LXBvb2wtcWMtUEhPNS5SbWRgKSwgSSBkaWQgUUMgYW5kIGV4cG9ydGVkIHRoZSBmaWx0ZXJlZCBkYXRhc2V0LiBIZXJlLCBJIHdpbGwgY29udGludWUgd29ya2luZyB3aXRoIHRoYXQgZGF0YXNldCBhbmQgYW5zd2VyIG91ciBiaW9sb2dpY2FsIHF1ZXN0aW9ucy4KCiMgR29hbAoKLSAgIEFuYWx5emUgdGhlIGZ1bGwgY2hpbWVyYSBzZXQgZmxvdyByZXN1bHRzIGZvciAqUEhPNXByKi1tQ2hlcnJ5IHJlcG9ydGVyLgotICAgRGV2ZWxvcCBhbiBhbmFseXNpcyBwaXBlbGluZSB0byBwZXJmb3JtIFFDLCBjb3JyZWN0aW9uIChpZiBuZWVkZWQpIGFuZCBwbG90dGluZyB0aGUgcmVzdWx0cy4KCiMgRGF0YQoKSW1wb3J0IHRoZSBiYWNrZ3JvdW5kIHN1YnRyYWN0ZWQgZGF0YQoKYGBge3J9CmRhdDAgPC0gcmVhZF90c3YoIi4uL2lucHV0LzIwMjMxMDIzLVBITzUtYmctc3VidHJhY3RlZC1kYXRhLnRzdiIsIGNvbF90eXBlcyA9ICJjY2NjZGRkZGRjIikKYGBgCgpGaWx0ZXIgdGhlIGRhdGEKCmBgYHtyfQpkYXQgPC0gZmlsdGVyKGRhdDAsIGhvc3QgIT0gIlBITzg0IiwgZmxhZyA9PSAicGFzcyIsIGRhdGUgIT0gIjAyLzEwIikgJT4lIAogICMgYmFzZWQgb24gcHJldmlvdXMgUUMsIHRoZSBmb2xsb3dpbmcgc2FtcGxlIChib3RoIHJlcGxpY2F0ZXMpIGhhdmUgaGlnaAogICMgdmFyaWFuY2UgLSBvbmUgYmlvbG9naWNhbCByZXBsaWNhdGUgaXMgaGlnaGx5IGV4cHJlc3NlZCwgd2hpbGUgdGhlIG90aGVyIAogICMgdHdvIGhhdmUgbU5lb24sIGJ1dCBiYXJlbHkgYW55IFJGUCBleHByZXNzaW9uLgogIG11dGF0ZSgKICAgIGhvc3QgPSBmY3RfcmVjb2RlKGhvc3QsIHBobzIgPSAicGhvMuKIhiIpLAogICAgZmxhZyA9IGlmZWxzZShwbGFzbWlkID09ICIyMzMiICYgaG9zdCA9PSAicGhvMiIsICJoaWdoLnZhciIsIGZsYWcpKQpgYGAKCk51bWJlciBvZiByZXBsaWNhdGVzIGxlZnQgZm9yIGVhY2ggc2FtcGxlCgpgYGB7cn0KZXhwdCA8LSBkYXQgJT4lIAogIGZpbHRlcihob3N0ICVpbiUgYygiUEhPMiIsICJwaG8yIiksICFwbGFzbWlkICVpbiUgYygiMTg4IiwgIjE5NCIpKSAlPiUgCiAgZ3JvdXBfYnkoZGF0ZSwgcGxhc21pZCwgaG9zdCkgJT4lIAogIHN1bW1hcml6ZShuID0gbigpLCAuZ3JvdXBzID0gImRyb3AiKQoKZXhwdCAlPiUgCiAgZ2dwbG90KGFlcyh4ID0gcGxhc21pZCwgeSA9IG4pKSArCiAgZ2VvbV9jb2woYWVzKGZpbGwgPSBob3N0KSkgKyAKICBmYWNldF9ncmlkKGRhdGUgfiAuKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygiUEhPMiIgPSAiZ3JheTMwIiwgInBobzIiID0gImdyYXk3MCIpKSArCiAgdGhlbWVfbWluaW1hbCgpICsgYmFja2dyb3VuZF9ncmlkKG1ham9yID0gIm5vbmUiKSArIHBhbmVsX2JvcmRlcihzaXplID0gMC41KSArCiAgc2NhbGVfeV9jb250aW51b3VzKG5hbWUgPSAiUmVwbGljYXRlcyIsIGJyZWFrcyA9IGMoNikpICsgeGxhYihOVUxMKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCksCiAgICAgICAgc3RyaXAudGV4dC55ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gMCksCiAgICAgICAgbGVnZW5kLnBvc2l0aW9uID0gInRvcCIpCmBgYAoKQ2hpbWVyYSBtYWtldXAgaW5mb3JtYXRpb24KCmBgYHtyfQptZXRhIDwtIHJlYWRfdHN2KCIuLi9pbnB1dC8yMDIzMDIwOC1jaGltZXJhLVBobzQtbWFrZXVwLnR4dCIsIGNvbF90eXBlcyA9ICJjY2NjYyIpCmBgYAoKIyMgU3VtbWFyaXplIGRhdGEKCkhlcmUgd2Ugd291bGQgbGlrZSBjYWxjdWxhdGUgdGhlIHJhdGlvIG9mIFJGUC9HRlAgZm9yIGVhY2ggY2hpbWVyYSAocGxhc21pZCkgYWNyb3NzIGFsbCByZXBsaWNhdGVzLCBpbmNsdWRpbmcgZnJvbSBkaWZmZXJlbnQgZGF5cy4gTm90ZSB0aGF0IHRoZSBwYXJhbWV0ZXIgb2YgaW50ZXJlc3QgaXMgYSByYXRpbywgd2hpY2ggY2FuIGJlIGVzdGltYXRlZCB1c2luZyBlaXRoZXIgIm1lYW5zIG9mIHJhdGlvcyIgb3IgInJhdGlvcyBvZiBtZWFucyIuIFRoZXNlIGFyZSBqdXN0IHR3byBzcGVjaWZpYyBpbnN0YW5jZXMgb2YgYSBtb3JlIGdlbmVyYWwgZXN0aW1hdG9yLCByZXByZXNlbnRpbmcgdHdvIGNob2ljZXMgb2YgdGhlIHdlaWdodHMuIFRoZSAibWVhbnMgb2YgcmF0aW9zIiBmaXJzdCBjYWxjdWxhdGVzIHRoZSByYXRpb3MgZm9yIGVhY2ggcmVwbGljYXRlIHdpdGhpbiBhIHBsYXNtaWQsIHRoZW4gYXZlcmFnZSB0aGVtLiBJbiB0aGlzIGNhbGN1bGF0aW9uLCBlYWNoIHJlcGxpY2F0ZSBpcyBnaXZlbiB0aGUgd2VpZ2h0IG9mIDEvbiAoZXF1YWwpLiBUaGUgInJhdGlvcyBvZiBtZWFucyIgZmlyc3Qgc3VtIHVwIHRoZSBHRlAgYW5kIFJGUCB2YWx1ZXMgc2VwYXJhdGVseSBhY3Jvc3MgdGhlIHJlcGxpY2F0ZXMgZm9yIGVhY2ggcGxhc21pZCwgdGhlbiB0YWtlIHRoZSByYXRpbyBiZXR3ZWVuIHRoZW0uIEluIHRoaXMgZXN0aW1hdG9yLCB0aGUgd2VpZ2h0IGZvciBlYWNoIHJlcGxpY2F0ZSBpcyB4IC8gc3VtKHgpLCB3aGVyZSB4IGlzIHRoZSBkZW5vbWluYXRvciBpbiB0aGUgcmF0aW8sIGkuZS4sIEdGUC4gSW4gb3RoZXIgd29yZHMsIHRoaXMgZXN0aW1hdG9yIHdpbGwgZ2l2ZSBtb3JlIHdlaWdodHMgdG8gdGhlIHJlcGxpY2F0ZXMgd2hlcmUgdGhlIGNoaW1lcmEgaGFkIGEgaGlnaGVyIGV4cHJlc3Npb24gbGV2ZWwuCgpCb3RoIGVzdGltYXRvcnMgYXJlIGtub3duIHRvIGJlIGJpYXNlZC4gV2Ugd2lsbCBpZ25vcmUgdGhhdCBmb3IgdGhlIG1vbWVudC4gSW4gdGVybXMgb2YgYSBjaG9pY2UgYmV0d2VlbiB0aGUgdHdvLCBpdCBzZWVtcyB0aGF0IHRoZXJlIGlzIG5vIHJlYXNvbiB0byBnaXZlIG1vcmUgd2VpZ2h0cyB0byB0aGUgZXhwZXJpbWVudHMgd2l0aCBhIGhpZ2hlciBHRlAgc2lnbmFsLiBTbywgdGhlICJtZWFucyBvZiByYXRpb3MiIHNlZW1zIGEgbW9yZSBuYXR1cmFsIGNob2ljZS4gSG93ZXZlciwgd2Ugd2lsbCBjYWxjdWx0YWUgYm90aCBhbmQgZGVjaWRlIGxhdGVyLgoKQSBmaW5hbCBxdWVzdGlvbiBpcyBob3cgdG8gY2FsY3VsYXRlIHRoZSB2YXJpYW5jZSBvZiB0aGUgcmF0aW8gZXN0aW1hdGUuIEFjY29yZGluZyB0byB0aGUgYHN1cnZleWAgcGFja2FnZSBbbWFudWFsXShodHRwczovL3JzdHVkaW8tcHVicy1zdGF0aWMuczMuYW1hem9uYXdzLmNvbS8xNzg5NjVfZmI2MGEwZjdiYmI0NGE2ZWEyMTk3MTNmYjFhODlhMjIuaHRtbCksIGFuIGFwcHJveGltYXRlIGVzdGltYXRvciBmb3IgdGhlIHZhcmlhbmNlIGlzCgokJApyID0gXGZyYWN7XGJhcnt5fX17XGJhcnt4fX0sIFx0ZXh0e3doZXJlfVwgXGJhcnt5fT1cZnJhY3sxfXtufVxzdW1fe2k9MX1ee259eV9pXCBcdGV4dHthbmR9XCBcYmFye3h9PVxmcmFjezF9e259XHN1bV97aT0xfV57bn14X2lcIFxcClxoYXR7Vn0ocikgPSAoMS1cZnJhY3tufXtOfSkoXGZyYWN7MX17XGJhcnt4fV4yfSlcZnJhY3tzX3JeMn17bn1cIFx0ZXh0e3doZXJlfVwgc19yXjI9XGZyYWN7MX17bi0xfVxzdW1fe2k9MX1ee259KHlfaS1yeF9pKV4yCiQkCgpBc3N1bWluZyB0aGF0IE5cPlw+biwgd2UgY2FuIGlnbm9yZSB0aGUgZmlyc3QgdGVybSBpbiB0aGUgdmFyaWFuY2UgZXN0aW1hdG9yLiBUaGUgcmVzdCBjYW4gYmUgY2FsY3VsYXRlZCBmcm9tIHRoZSBkYXRhCgpgYGB7cn0KZGF0c3VtIDwtIGRhdCAlPiUKICBmaWx0ZXIoIWlzLm5hKHBsYXNtaWQpKSAlPiUgCiAgZ3JvdXBfYnkocGxhc21pZCwgaG9zdCkgJT4lIAogIHN1bW1hcml6ZSgKICAgICBuID0gbigpLAogICAgbUcgPSBtZWFuKEJMMS5IKSwKICAgIG1SID0gbWVhbihZTDIuSCksCiAgICAgQSA9IG1lYW4oWUwyLkgvQkwxLkgpLAogICAgIHIgPSBtUi9tRywKICAgIHMyID0gMS8obi0xKSpzdW0oKFlMMi5IIC0gcipCTDEuSCleMiksCiAgICB2ciA9IDEvKG1HXjIpKnMyL24sCiAgICBzZSA9IHNxcnQodnIpLAogICAgLmdyb3VwcyA9ICJkcm9wIgogICkgJT4lIAogIHNlbGVjdCgtczIsIC1yLCAtdnIpIyAlPiUgCiAgI3Bpdm90X3dpZGVyKG5hbWVzX2Zyb20gPSBob3N0LCB2YWx1ZXNfZnJvbSA9IEJMMS5IOmBuUi9HYCkgJT4lIAogICNtdXRhdGUoYHBobzLiiIYvUEhPMmAgPSBgUi9HX3BobzLiiIZgL2BSL0dfUEhPMmAsCiAgIyAgICAgICBgbi5waG8y4oiGL1BITzJgID0gYG5SL0dfcGhvMuKIhmAvYG5SL0dfUEhPMmApCmBgYAoKRm9yIGVhY2ggY2hpbWVyYSwgd2Ugd291bGQgYWxzbyBsaWtlIHRvIGNhbGN1bGF0ZSAqKnRocmVlIHZhbHVlcyoqOgoKMS4gIEEgaW4gKnBobzLiiIYqOiB0aGlzIGlzIGl0cyBiYXNlIGFjdGl2aXR5IHdpdGhvdXQgUGhvMgoyLiAgQSBpbiAqUEhPMio6IHRoaXMgaXMgaXRzIGZ1bGwgYWN0aXZpdHkgd2l0aCBQaG8yCjMuICBBX1BITzIgLyBBX3BobzLiiIY6IHRoaXMgaXMgdGhlIFBobzIgZW5oYW5jZW1lbnQgb2YgYWN0aXZpdHkKCldlIGFzc2lnbiB0aGUgY2hpbWVyYXMgaW50byBzZXZlcmFsIGdyb3VwcywgYmFzZWQgb24gdGhlaXIgQV9QSE8yIGFuZCBBX1BITzIvQV9waG8y4oiGCgpgYGB7cn0KeGltZXJhIDwtIGRhdHN1bSAlPiUKICBwaXZvdF93aWRlcihpZF9jb2xzID0gcGxhc21pZCwgbmFtZXNfZnJvbSA9IGhvc3QsCiAgICAgICAgICAgICAgdmFsdWVzX2Zyb20gPSBjKEEsIHNlKSkgJT4lIAogIG11dGF0ZSgKICAgIHJBX1BITzIgPSBBX1BITzIgLyBBX1BITzJbcGxhc21pZCA9PSAiMTk0Il0sCiAgICByQV9waG8yID0gQV9waG8yIC8gQV9waG8yW3BsYXNtaWQgPT0gIjE5NCJdLAogICAgYm9vc3QgPSBBX1BITzIgLyBBX3BobzIsCiAgICBncm91cCA9IGNhc2Vfd2hlbigKICAgICAgcGxhc21pZCAlaW4lIGMoIjE4OCIsICIxOTQiKSB+ICJyZWYiLAogICAgICByQV9QSE8yIDwgMC4yICAgICAgICAgICAgICAgIH4gIm4uZi4iLAogICAgICAuZGVmYXVsdCA9ICJjaGltZXJhIgogICAgKSwKICAgIGdyb3VwID0gZmN0X3JlbGV2ZWwoZ3JvdXAsICJyZWYiLCAiY2hpbWVyYSIsICJuLmYuIikKICApICU+JSAKICBsZWZ0X2pvaW4oc2VsZWN0KG1ldGEsIHBsYXNtaWQsIHNldCwgc3ltYm9sKSwgYnkgPSAicGxhc21pZCIpICU+JSAKICBtdXRhdGUoc3ltYm9sID0gZmN0X3Jlb3JkZXIoc3ltYm9sLCByQV9QSE8yLCAuZGVzYyA9IFRSVUUpKSAlPiUgCiAgcmVsb2NhdGUoYyhzZXQsIHN5bWJvbCwgZ3JvdXApLCAuYWZ0ZXIgPSBwbGFzbWlkKQpgYGAKCkV4cG9ydCB0aGUgc3VtbWFyaXplZCBkYXRhCgpgYGB7cn0Kd3JpdGVfdHN2KHhpbWVyYSwgZmlsZSA9ICIuLi9vdXRwdXQvMjAyMzExMjUtUEhPNXByLWNoaW1lcmEtc3VtbWFyaXplZC50c3YiKQpgYGAKClRvIGJlIGFibGUgdG8gcGxvdCBhbGwgdGhlIGRhdGEgcG9pbnRzLCBsZXQncyBnZW5lcmF0ZSBhbm90aGVyIGRhdGEgZnJhbWUgd2l0aCB0aGUgaW5kaXZpZHVhbCByYXRpb3MuCgpgYGB7cn0KZGF0X3NlcCA8LSBkYXQgJT4lCiAgZmlsdGVyKCFpcy5uYShwbGFzbWlkKSkgJT4lIAogIG11dGF0ZShBID0gWUwyLkgvQkwxLkgpICU+JSAKICBzZWxlY3QocGxhc21pZCwgaG9zdCwgQkwxLkgsIFlMMi5ILCBBLCBmbGFnKSAlPiUgCiAgbGVmdF9qb2luKHNlbGVjdChtZXRhLCBwbGFzbWlkLCBzZXQsIHN5bWJvbCksIGJ5ID0gInBsYXNtaWQiKSMgJT4lIAogICNtdXRhdGUoc3ltYm9sID0gZmN0X3Jlb3JkZXIoc3ltYm9sLCByQV9QSE8yLCAuZGVzYyA9IFRSVUUpKSMgJT4lIApgYGAKCiMgQW5hbHlzaXMKCiMjIFBsb3R0aW5nIGZ1bmN0aW9ucwoKPCEtLS0KU2V0IHVwIGNvbW1vbiBwYXJhbWV0ZXJzIGZvciB0aHJlc2hvbGRpbmcgYW5kIHBsb3R0aW5nCmBgYHtyfQojIHJlZmVyZW5jZSBQaG80IHBsYXNtaWQgaWRzCnJlZnMgPC0gYygiMTg4IiwgIjE5NCIpCiMgY29sb3JzCmRhdGUuY29sb3JzID0gYyhicmV3ZXIucGFsKG5hbWU9IkRhcmsyIiwgbiA9IDgpLCBicmV3ZXIucGFsKG5hbWU9IlBhaXJlZCIsIG4gPSA4KSkKaG9zdC5jb2xvcnMgPSBjKCJQSE8yIiA9ICJncmF5MzAiLCAicGhvMiIgPSAiZ3JheTcwIikKcG9pbnQuY29sb3JzID0gYygiUEhPMiIgPSAiZm9yZXN0Z3JlZW4iLCAicGhvMiIgPSAicHVycGxlNCIpCiMgCmBgYAotLS0+CgpTb3VyY2UgdGhlIHNjcmlwdHMKCmBgYHtyfQpzb3VyY2UoIi4uL3NjcmlwdC8yMDI0MDIxMS1jaGltZXJhLWRhdGEtcGxvdHRpbmctZnVuY3Rpb25zLlIiKQpgYGAKCk1vZGlmeSB0aGUgY29tcG9uZW50IHBsb3R0aW5nIGZ1bmN0aW9uIGZvciBzcGVjaWFsIHB1cnBvc2VzCgpgYGB7cn0KbXlfcGxvdF9yYXRpbyA8LSBmdW5jdGlvbihzZWxlY3Rpb24pewogICMgY3VzdG9tIGNvbG9ycyBmb3IgdGhpcyBmdW5jdGlvbgogIGRhdGUuY29sb3JzID0gYyhicmV3ZXIucGFsKG5hbWU9IkRhcmsyIiwgbiA9IDgpLCBicmV3ZXIucGFsKG5hbWU9IlBhaXJlZCIsIG4gPSA4KSkKICBob3N0LmNvbG9ycyA9IGMoIlBITzIiID0gImdyYXkzMCIsICJwaG8yIiA9ICJncmF5NzAiKQogIHBvaW50LmNvbG9ycyA9IGMoIlBITzIiID0gImZvcmVzdGdyZWVuIiwgInBobzIiID0gInB1cnBsZTQiKQogICMgcHJlcGFyZSBkYXRhCiAgdG1wIDwtIG15X2RhdGFfcHJlcChzZWxlY3Rpb24pCiAgIyBwbG90dGluZwogIHAgPC0gdG1wICU+JSAKICAgIHNlbGVjdCgtYyhGU0MuSCwgbkdGUCwgblJGUCwgZmxhZykpICU+JSAKICAgIG11dGF0ZShgUi9HYCA9IFlMMi5IL0JMMS5IKSAlPiUgCiAgICBwaXZvdF9sb25nZXIoY29scyA9IGMoQkwxLkgsIFlMMi5ILCBgUi9HYCksIAogICAgICAgICAgICAgICAgIG5hbWVzX3RvID0gInBhcmFtZXRlciIsIHZhbHVlc190byA9ICJ2YWx1ZSIpICU+JSAKICAgIG11dGF0ZShwYXJhbWV0ZXIgPSBmYWN0b3IocGFyYW1ldGVyLCBsZXZlbHMgPSBjKCJSL0ciLCAiWUwyLkgiLCAiQkwxLkgiKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbGFiZWxzID0gYygiUkZQL0dGUCIsICJQSE81cFJGUCIsICJQaG80LUdGUCIpKSkgJT4lIAogICAgZ2dwbG90KGFlcyh4ID0gc3ltYm9sLCB5ID0gdmFsdWUsIGdyb3VwID0gaG9zdCkpICsgCiAgICBzdGF0X3N1bW1hcnkoYWVzKGdyb3VwID0gaG9zdCksIGZ1bi5kYXRhID0gIm1lYW5fY2xfYm9vdCIsIGdlb20gPSAiZXJyb3JiYXIiLAogICAgICAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2UoMC41KSwgd2lkdGggPSAwLjMpICsKICAgIGdlb21fYmFyKGFlcyhmaWxsID0gaG9zdCksIHdpZHRoID0gMC41LCBhbHBoYSA9IDAuOCwKICAgICAgICAgICAgIHN0YXQgPSAic3VtbWFyeSIsIGZ1biA9ICJtZWFuIiwgcG9zaXRpb24gPSBwb3NpdGlvbl9kb2RnZSgwLjUpKSArCiAgICBnZW9tX3BvaW50KGRhdGEgPSBmdW5jdGlvbih4KSBzdWJzZXQoeCwgIXN5bWJvbCAlaW4lIGMoIkNDQ0NDIiwgIlNTU1NTIikpLAogICAgICAgICAgICAgICBhZXMoZ3JvdXAgPSBob3N0LCBjb2xvciA9IGRhdGUpLCBzaXplID0gMSwgc2hhcGUgPSAzLCBhbHBoYSA9IDAuOSwKICAgICAgICAgICAgICAgcG9zaXRpb24gPSBwb3NpdGlvbl9qaXR0ZXJkb2RnZShkb2RnZS53aWR0aCA9IDAuNSwgaml0dGVyLndpZHRoID0gMC4xKSkgKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGRhdGUuY29sb3JzLCBndWlkZSA9ICJub25lIikgKwogICAgI2dlb21fcG9pbnQoZGF0YSA9IGZ1bmN0aW9uKHgpIHN1YnNldCh4LCAhc3ltYm9sICVpbiUgYygiQ0NDQ0MiLCAiU1NTU1MiKSksCiAgICAjICAgICAgICAgICBhZXMoZ3JvdXAgPSBob3N0LCBjb2xvciA9IGhvc3QpLCBzaXplID0gMSwgc2hhcGUgPSAzLCBhbHBoYSA9IDAuOSwKICAgICMgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25faml0dGVyZG9kZ2UoZG9kZ2Uud2lkdGggPSAwLjUsIGppdHRlci53aWR0aCA9IDAuMSkpICsKICAgICNzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gcG9pbnQuY29sb3JzKSArCiAgICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBob3N0LmNvbG9ycykgKwogICAgZmFjZXRfZ3JpZChwYXJhbWV0ZXJ+Z3JvdXAsIHNjYWxlcyA9ICJmcmVlIiwgc3BhY2UgPSAiZnJlZV94IikgKwogICAgdGhlbWVfYncoYmFzZV9zaXplID0gMTgpICsgYmFja2dyb3VuZF9ncmlkKG1pbm9yID0gIm5vbmUiKSArIAogICAgeGxhYigiUGhvNCBjaGltZXJhIikgKwogICAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSAzMCwgaGp1c3QgPSAxLCBmYW1pbHkgPSAibW9ubyIpLAogICAgICAgICAgbGVnZW5kLnBvc2l0aW9uID0gInRvcCIsCiAgICAgICAgICBheGlzLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKQogIHJldHVybihwKSAKfQpgYGAKCiMjIFBobzQgY2hpbWVyYSBwcm90ZWluIGxldmVsIHZhcmlhdGlvbgoKV2hhdCBpcyB0aGUgZGlzdHJpYnV0aW9uIG9mIFBobzQgY2hpbWVyYSBwcm90ZWluIGxldmVscz8gSXMgdGhlIG1DaGVycnkvbU5lb24gcmF0aW8gYSBmYWl0aGZ1bCBtZWFzdXJlIG9mIHRoZSBjaGltZXJhJ3MgYWN0aXZpdGllcz8KCkRpc3RyaWJ1dGlvbiBvZiBQaG80LW1OZW9uIGxldmVscyBncm91cGVkIGJ5IHBsYXNtaWQgYW5kIGhvc3QuCgpgYGB7cn0KaG9zdC5sYWJlbHMgPSBjKCJQSE8yIiwgInBobzLiiIYiKQpwb2ludC5jb2xvcnMgPSBjKCJQSE8yIiA9ICJmb3Jlc3RncmVlbiIsICJwaG8yIiA9ICJwdXJwbGU0IikKcDEgPC0gZGF0ICU+JSAKICBmaWx0ZXIoIWlzLm5hKHBsYXNtaWQpKSAlPiUgCiAgbXV0YXRlKHBsYXNtaWQgPSBmY3RfcmVvcmRlcihwbGFzbWlkLCBCTDEuSCwgLmZ1biA9IG1lZGlhbikgJT4lIAogICAgICAgICAgIGZjdF9yZWxldmVsKCIxOTQiLCAiMTg4IikpICU+JSAKICBnZ3Bsb3QoYWVzKHggPSBwbGFzbWlkLCB5ID0gQkwxLkgpKSArCiAgZ2VvbV9wb2ludChhZXMoY29sb3IgPSBob3N0KSwgcG9zaXRpb24gPSBwb3NpdGlvbl9qaXR0ZXIoMC4xKSwKICAgICAgICAgICAgIHNpemUgPSAxLjEpICsgCiAgc2NhbGVfY29sb3JfbWFudWFsKCJIb3N0IiwgdmFsdWVzID0gcG9pbnQuY29sb3JzLCBsYWJlbHMgPSBob3N0LmxhYmVscykgKwogIHNjYWxlX3lfbG9nMTAoYnJlYWtzID0gYygxMDAsIDEwMDAsIDEwMDAwKSwgZXhwYW5kID0gZXhwYW5zaW9uKG11bHQgPSAwLjEpKSArCiAgc2NhbGVfeF9kaXNjcmV0ZShleHBhbmQgPSBleHBhbnNpb24obXVsdCA9IDAuMDMpKSArCiAgeGxhYigiUGhvNCBjb25zdHJ1Y3RzIikgKyB5bGFiKCJQaG80LW1OZW9uIChhLnUuKSIpICsKICB0aGVtZV9jb3dwbG90KCkgKyBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIiwgc2l6ZSA9IDEuMikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHNpemUgPSByZWwoMC42KSwgdmp1c3QgPSAwLjUpLAogICAgICAgIGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMC44KSksCiAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDAuOSkpLAogICAgICAgIGF4aXMubGluZSA9IGVsZW1lbnRfYmxhbmsoKSwKICAgICAgICBsZWdlbmQucG9zaXRpb24gPSBjKDAuMDUsIDAuOSksCiAgICAgICAgbGVnZW5kLmRpcmVjdGlvbiA9ICJob3Jpem9udGFsIiwKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChmYWNlID0gMykpCnAxCiNnZ3NhdmUoIi4uL2ltZy8yMDI0MDMwNy1QaG80LWNoaW1lcmEtcHJvdGVpbi1sZXZlbC12YXJpYXRpb24ucG5nIiwgCiMgICAgICAgd2lkdGggPSA2LCBoZWlnaHQgPSAzKQpgYGAKCkFsbCBjaGltZXJhIChTY1BobzQgYW5kIENnUGhvNCByZW1vdmVkKSBwcm90ZWluIGxldmVscyBieSBob3N0CgpgYGB7cn0KdG1wIDwtIGRhdCAlPiUgCiAgZmlsdGVyKHBsYXNtaWQgPT0gIjE5NCIsIGhvc3QgPT0gIlBITzIiKQoKbG0gPC0gbG0oWUwyLkggfiBCTDEuSCwgZGF0YSA9IHRtcCkKc3VtbWFyeShsbSkKCnAyIDwtIHRtcCAlPiUgCiAgZ2dwbG90KGFlcyh4ID0gQkwxLkgsIHkgPSBZTDIuSCkpICsKICBnZW9tX3BvaW50KHNpemUgPSAxLjUpICsgCiAgc3RhdF9zbW9vdGgobWV0aG9kID0gImxtIiwgZm9ybXVsYSA9IHkgfiB4KSArCiAgeGxhYigiUGhvNC1tTmVvbiIpICsgeWxhYigiUEhPNXByLW1DaGVycnkiKSArCiAgdGhlbWVfY293cGxvdCgpICsgCiAgcGFuZWxfYm9yZGVyKGNvbG9yID0gImdyYXkzMCIsIHNpemUgPSAxLjIpICsKICB0aGVtZShheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMikpCiAgKQoKcDIKZ2dzYXZlKCIuLi9pbWcvMjAyNDAzMDctU2NQaG80LW1DaGVycnktdnMtbU5lb24tY29uc2lzdGVudC5wbmciLCAKICAgICAgIHdpZHRoID0gNCwgaGVpZ2h0ID0gMy41KQpgYGAKCmBgYHtyfQpob3N0LmNvbG9ycyA9ICBjKCJQSE8yIiA9ICJncmF5NjAiLCAicGhvMiIgPSAib3JhbmdlIikKCnRtcCA8LSBkYXQgJT4lIAogIGZpbHRlcihwbGFzbWlkICVpbiUgYygiMTg4IiwgIjE5NCIpLCAKICAgICAgICAgZGF0ZSAlaW4lIGMoIjAyLzA4IiwgIjAyLzExIiwgIjAyLzE4IiwgIjAyLzIxIiwgIjAyLzIzIiwgIjAzLzMxIikpICU+JSAKICBtdXRhdGUoQSA9IFlMMi5IL0JMMS5ILAogICAgICAgICBQaG80ID0gZmFjdG9yKHBsYXNtaWQsIGxldmVscyA9IGMoIjE5NCIsICIxODgiKSwgCiAgICAgICAgICAgICAgICAgICAgICAgbGFiZWxzID0gYygiU2NQaG80IiwgIkNnUGhvNCIpKSkgCnAzIDwtIHRtcCAlPiUgCiAgZ2dwbG90KGFlcyh4ID0gZGF0ZSwgeSA9IEEpKSArIAogIGdlb21fYmFyKGFlcyhmaWxsID0gaG9zdCksIHN0YXQgPSAic3VtbWFyeSIsIGZ1biA9ICJtZWFuIiwgCiAgICAgICAgICAgcG9zaXRpb24gPSBwb3NpdGlvbl9kb2RnZSgwLjkpLCBhbHBoYSA9IDAuOSkgKwogIGdlb21fcG9pbnQoYWVzKGdyb3VwID0gaG9zdCksIHNpemUgPSAwLjYsIHNoYXBlID0gMywKICAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpKSArIAogIHNjYWxlX2ZpbGxfbWFudWFsKCJIb3N0IiwgdmFsdWVzID0gaG9zdC5jb2xvcnMsIGxhYmVscyA9IGhvc3QubGFiZWxzKSArCiAgc2NhbGVfeF9kaXNjcmV0ZShsYWJlbHMgPSAxOjYpICsKICAjc3RhdF9zdW1tYXJ5KGZ1bi5kYXRhID0gIm1lYW5fc2UiLCBnZW9tID0gInBvaW50cmFuZ2UiLCBjb2xvciA9ICJyZWQiKSArCiAgZmFjZXRfZ3JpZChQaG80IH4gLikgKwogIHlsYWIoIm1DaGVycnkvbU5lb24iKSArIHhsYWIoIlJlcGxpY2F0ZSIpICsKICB0aGVtZV9jb3dwbG90KCkgKyBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIiwgc2l6ZSA9IDEuMikrCiAgdGhlbWUoYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMC43KSksCiAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEpKSwKICAgICAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgc3RyaXAuYmFja2dyb3VuZCA9IGVsZW1lbnRfYmxhbmsoKSwKICAgICAgICBsZWdlbmQucG9zaXRpb24gPSAidG9wIiwKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgwLjkpKSwKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDAuOCksIGZhY2UgPSAzKSkKcDMKI2dnc2F2ZSgiLi4vaW1nLzIwMjQwMzA3LUNnUGhvNC1tQ2hlcnJ5LXZzLW1OZW9uLWNvbnNpc3RlbnQucG5nIiwgCiMgICAgICAgd2lkdGggPSA0LCBoZWlnaHQgPSAzLjIpCgojIHNhbXBsZSBzaXplIHBlciBkYXkgb2YgZXhwZXJpbWVudAp0bXAgJT4lIGNvdW50KGRhdGUsIFBobzQsIGhvc3QpICU+JSAKICBwaXZvdF93aWRlcihuYW1lc19mcm9tID0gaG9zdCwgdmFsdWVzX2Zyb20gPSBuKQpgYGAKCiMjIEhpZ2ggdmFyaWFuY2Ugc2FtcGxlcwoKU3VtbWFyaXplIHRoZSBiYWNrZ3JvdW5kIHN1YnRyYWN0ZWQgZGF0YSBieSBjYWxjdWxhdGluZyB0aGUgbWVhbnMgYW5kIGN2IGZvciBlYWNoIHN0cmFpbi4KCmBgYHtyfQpjdiA8LSBkYXQgJT4lIAogIHNlbGVjdCgtbkdGUCwgLW5SRlApICU+JQogIHBpdm90X2xvbmdlcihGU0MuSDpZTDIuSCwgbmFtZXNfdG8gPSAicGFyYW1ldGVyIiwgdmFsdWVzX3RvID0gImludGVuc2l0eSIpICU+JSAKICBncm91cF9ieShkYXRlLCBwbGFzbWlkLCBob3N0LCBwYXJhbWV0ZXIpICU+JSAKICBzdW1tYXJpemUoCiAgICBuID0gbigpLAogICAgbWVhbiA9IG1lYW4oaW50ZW5zaXR5KSwKICAgIGN2ID0gc2QoaW50ZW5zaXR5KS9tZWFuKGludGVuc2l0eSksCiAgICAuZ3JvdXBzID0gImRyb3AiCiAgKSAlPiUgCiAgYXJyYW5nZShkZXNjKGN2KSkKYGBgCgpVc2UgdGhlIGNvbnRyb2wgc3RyYWluIChwSDE5NCB3aXRoIFBITzIpIHRvIGlkZW50aWZ5IGFuZCBjb3JyZWN0IGZvciBzeXN0ZW1hdGljIGJpYXNlcwoKYGBge3J9CmNvbnRyb2wgPC0gZmlsdGVyKGRhdCwgcGxhc21pZCA9PSAiMTk0IiwgaG9zdCA9PSAiUEhPMiIpICU+JSAKICBzZXBhcmF0ZSh3ZWxsLCBpbnRvID0gYygicm93IiwgImNvbCIpLCBzZXAgPSAxKSAlPiUgCiAgZHJvcGxldmVscygpCmBgYAoKTW9kZWwgZm9yIG1OZW9uCgpgYGB7cn0KZ2ZwLm1vZGVsLjAgPC0gbG0oQkwxLkggfiBsb2cxMChldmVudHMpICsgZGF0ZSArIHJvdypjb2wsIGRhdGEgPSBjb250cm9sKQpzdGVwKGdmcC5tb2RlbC4wKQpnZnAubW9kZWwuMSA8LSBsbShCTDEuSCB+IGRhdGUgKyBjb2wsIGRhdGEgPSBjb250cm9sKQpgYGAKCk1vZGVsIGZvciBQSE81cHI6OlJGUAoKYGBge3J9CnJmcC5tb2RlbC4wIDwtIGxtKFlMMi5IIH4gbG9nMTAoZXZlbnRzKSArIGRhdGUgKyByb3cqY29sLCBkYXRhID0gY29udHJvbCkKc3RlcChyZnAubW9kZWwuMCkKcmZwLm1vZGVsLjEgPC0gbG0oWUwyLkggfiBsb2cxMChldmVudHMpICsgZGF0ZSArIHJvdyArIGNvbCwgZGF0YSA9IGNvbnRyb2wpCmBgYAoKPiB0aGVyZSBhcmUgbW9yZSBzeXN0ZW1hdGljIHNoaWZ0cyBpbiB0aGUgUkZQLCBzaWduaWZpY2FudCBmb3Igcm93LCBjb2wsIGRhdGUgYW5kIGFsc28gXCMgb2YgZXZlbnRzIGhvd2V2ZXIsIEkgd29uJ3QgYmUgcmVtb3ZpbmcgdGhlc2UgZWZmZWN0cyB5ZXQsIGJlY2F1c2UgSSd2ZSBmb3VuZCB0aGF0IFJGUC9HRlAgcmF0aW9zIGFyZSBwcmV0dHkgY29uc2lzdGVudCBhY3Jvc3MgZGF5cy4gSW4gb3RoZXIgd29yZHMsIHRoZSB2YXJpYXRpb24gaW4gR0ZQIGFuZCBSRlAgbWF5IGJlIGNhbmNlbGxlZCBvdXQuCgpDaGVjayBmb3IgZWFjaCBwbGFzbWlkIGhvdyBjb25zaXN0ZW50IGFyZSB0aGUgbWVhc3VyZW1lbnRzIGJldHdlZW4gZGF5cwoKYGBge3J9CnRtcCA8LSBkYXQgJT4lIAogICMgcmVtb3ZlIG9uZSBzYW1wbGUgd2l0aCBvbmx5IG9uZSB2YWxpZCBkYXkgb2YgZXhwZXJpbWVudAogIGZpbHRlcighKHBsYXNtaWQgPT0gIjIxOCIgJiBob3N0ID09ICJQSE8yIiksICFwbGFzbWlkICVpbiUgYygiMTg4IiwgIjE5NCIsIE5BKSkgJT4lIAogIG5lc3QoZGF0YSA9IGMoZGF0ZSwgQkwxLkgsIFlMMi5IKSwgLmJ5ID0gYyhwbGFzbWlkLCBob3N0KSkKCmRheS52YXIuZ2ZwIDwtIHRtcCAlPiUgCiAgbXV0YXRlKG1vZGVsID0gbWFwKGRhdGEsIGZ1bmN0aW9uKGRmKSBsbShCTDEuSCB+IGRhdGUsIGRhdGEgPSBkZikpLAogICAgICAgICB0aWRpZWQgPSBtYXAobW9kZWwsIGJyb29tOjp0aWR5KSkgJT4lIAogIHVubmVzdCh0aWRpZWQpICU+JSAKICBmaWx0ZXIodGVybSAhPSAiKEludGVyY2VwdCkiKSAlPiUgCiAgbXV0YXRlKHAuYWRqID0gcC5hZGp1c3QocC52YWx1ZSwgbWV0aG9kID0gIkJIIikpICU+JSAKICBzZWxlY3QoLWRhdGEsIC1tb2RlbCkgJT4lIAogIGZpbHRlcihwLmFkaiA8IDAuMTApICU+JSAKICBhcnJhbmdlKHBsYXNtaWQsIGhvc3QpCgpkYXkudmFyLnJmcCA8LSB0bXAgJT4lIAogIG11dGF0ZShtb2RlbCA9IG1hcChkYXRhLCBmdW5jdGlvbihkZikgbG0oWUwyLkggfiBkYXRlLCBkYXRhID0gZGYpKSwKICAgICAgICAgdGlkaWVkID0gbWFwKG1vZGVsLCBicm9vbTo6dGlkeSkpICU+JSAKICB1bm5lc3QodGlkaWVkKSAlPiUgCiAgZmlsdGVyKHRlcm0gIT0gIihJbnRlcmNlcHQpIikgJT4lIAogIG11dGF0ZShwLmFkaiA9IHAuYWRqdXN0KHAudmFsdWUsIG1ldGhvZCA9ICJCSCIpKSAlPiUgCiAgc2VsZWN0KC1kYXRhLCAtbW9kZWwpICU+JSAKICBmaWx0ZXIocC5hZGogPCAwLjEwKSAlPiUgCiAgYXJyYW5nZShwbGFzbWlkLCBob3N0KQpgYGAKCmBgYHtyfQojIGV4dHJhY3QgeGltZXJhIG5hbWVzCnJlZnMgPC0gYygiMTg4IiwiMTk0IikKIyBtYWtlIGEgdGVzdCBzZXQKZGF5LnZhci5nZnAubGlzdCA8LSB1bmlxdWUoZGF5LnZhci5nZnAkcGxhc21pZCkKZGF5LnZhci5yZnAubGlzdCA8LSB1bmlxdWUoZGF5LnZhci5yZnAkcGxhc21pZCkKYGBgCgpIaWdoIGRheS10by1kYXkgR0ZQIHZhcmlhbmNlOiBgciBkYXkudmFyLmdmcC5saXN0YCBIaWdoIGRheS10by1kYXkgUkZQIHZhcmlhbmNlOiBgciBkYXkudmFyLnJmcC5saXN0YAoKUGxvdHRpbmcgY29tcG9uZW50cyBmb3IgY2hpbWVyYXMgd2l0aCBoaWdoIGRheS10by1kYXkgdmFyaWFuY2UgaW4gUGhvNC1tTmVvbgoKYGBge3J9CnAgPC0gbXlfcGxvdF9yYXRpbyhjKHJlZnMsZGF5LnZhci5nZnAubGlzdCkpIyArIApwCmBgYAoKPiBXYXRjaCBvdXQgZm9yIENTQ3NjQywgU0NDc1MsIFNDQ3NTCgpQbG90dGluZyBjb21wb25lbnRzIGZvciBjaGltZXJhcyB3aXRoIGhpZ2ggZGF5LXRvLWRheSB2YXJpYW5jZSBpbiAqUEhPNXByKi1tQ2hlcnJ5CgpgYGB7cn0KcCA8LSBteV9wbG90X3JhdGlvKGMocmVmcyxkYXkudmFyLnJmcC5saXN0KSkKcApgYGAKCj4gbW9zdCBvZiB0aGUgZGF5LXRvLWRheSB2YXJpYW5jZSBhcmUgY2FuY2VsZWQgb3V0IGFmdGVyIFJGUC9HRlAgbm9ybWFsaXphdGlvbgoKIyMgQWxsIGNoaW1lcmEsIHNjYXR0ZXIgcGxvdAoKUGxvdCBhbGwgY2hpbWVyYXMsIGNvbG9yaW5nIGJhc2VkIG9uIFAySUQgc291cmNlCmBgYHtyfQpteV9zY2F0dGVyX3Bsb3RfZml4IDwtIGZ1bmN0aW9uKCl7CiAgIyB0aGlzIGZ1bmN0aW9uIGlzIHRoZSBzYW1lIGFzIHRoZSBvbmUgaW4gdGhlIHNjcmlwdCBmaWxlLCBidXQgaXMgdXNlZCB0byAKICAjIHBsb3QgcmVnaW9uIDQgZWZmZWN0cyBhbG9uZSwgYW5kIGRvZXNuJ3QgdGFrZSBhbnkgaW5wdXQKICBzMSA9IG15X2RhdGFfc2VsZWN0KHBhdHRlcm4gPSAiWFhYQ0NYIikKICBzMiA9IG15X2RhdGFfc2VsZWN0KHBhdHRlcm4gPSAiWFhYU1NYIikKICBzY2F0dGVyLmNvbG9ycyA9IGMoIlNjUGhvNCIgPSAiZm9yZXN0Z3JlZW4iLCAiQ2dQaG80IiA9ICJibHVlMyIsIAogICAgICAgICAgICAgICAgICAgICAiUDJJRDpDZyIgPSAiZGVlcHNreWJsdWUiLCAiUDJJRDpTYyIgPSAicGFsZWdyZWVuMiIsCiAgICAgICAgICAgICAgICAgICAgICJQMklEOm1peGVkIiA9ICJncmF5MjAiKQogIHNjYXR0ZXIuc2l6ZSA9IGMoIlNjUGhvNCIgPSAzLjUsICJDZ1BobzQiID0gMy41LAogICAgICAgICAgICAgICAgICAgICAiUDJJRDpDZyIgPSAyLjUsICJQMklEOlNjIiA9IDIuNSwgIlAySUQ6bWl4ZWQiID0gMi41KQogIHAgPC0geGltZXJhICU+JSAKICAgICMgZXhjbHVkZSB0aGUgYWx0ZXJuYXRpdmUgYnJlYWsgcG9pbnQgc2V0cyAiQSIgYW5kICJCIgogICAgIyBpbiBwYXJ0aWN1bGFyLCBwSDI5NCBpcyBhbiBhbHRlcm5hdGl2ZSBicmVhayBwb2ludCBDQ0Njcywgd2hlcmUgUDJJRDpDZwogICAgIyBleHRlbmRzIHRvIGFhIDI3MCBpbnN0ZWFkIG9mIDQ1OC4gaXQgaGFzIG5lYXJseSB0aGUgc2FtZSBhY3Rpdml0aWVzIGFzCiAgICAjIENnUGhvNC4gSG93ZXZlciwgdGhlIENDQ0NTIGluIHRoZSBtYWluIHNldCBoYXMgc2lnbmlmaWNhbnRseSByZWR1Y2VkCiAgICAjIGFjdGl2aXRpZXMgYm90aCB3aXRoIGFuZCB3aXRob3V0IFBobzQuIFdlIGxhdGVyIHRlc3RlZCB3aGV0aGVyIHRoZSAKICAgICMgYWRkaXRpb25hbCBQMklEOkNnNCByZXNjdWVzIHRoZSBlZmZlY3QgKHNlZSBDZzRleHQgYmVsb3cpIGFuZCBpdCBkaWRuJ3QKICAgICMgc28gZmFyLCB0aGlzIHNlZW1zIHRvIGJlIGFuIG9uZS1vZmYuIHdlIG5lZWQgdG8gZnVydGhlciBpbnZlc3RpZ2F0ZSBpdHMKICAgICMgYWN0aXZpdGllcy4KICAgIGZpbHRlcihzZXQgJWluJSBjKCJNIiwgIlMiKSkgJT4lIAogICAgbXV0YXRlKEFfUEhPMiA9IHNpZ25pZihBX1BITzIsIGRpZ2l0cyA9IDIpLAogICAgICAgICAgIEFfcGhvMiA9IHNpZ25pZihBX3BobzIsIGRpZ2l0cyA9IDIpLAogICAgICAgICAgIGdyb3VwID0gY2FzZV93aGVuKAogICAgICAgICAgICAgc3ltYm9sID09ICJDQ0NDQyIgfiAiQ2dQaG80IiwKICAgICAgICAgICAgIHN5bWJvbCA9PSAiU1NTU1MiIH4gIlNjUGhvNCIsCiAgICAgICAgICAgICBwbGFzbWlkICVpbiUgczEgfiAiUDJJRDpDZyIsCiAgICAgICAgICAgICBwbGFzbWlkICVpbiUgczIgfiAiUDJJRDpTYyIsCiAgICAgICAgICAgICAuZGVmYXVsdCA9ICJQMklEOm1peGVkIgogICAgICAgICAgICksCiAgICAgICAgICAgZ3JvdXAgPSBmY3RfcmVsZXZlbChncm91cCwgbmFtZXMoc2NhdHRlci5jb2xvcnMpKSkgJT4lIAogICAgZ2dwbG90KGFlcyh4ID0gQV9QSE8yLCB5ID0gQV9waG8yLCBsYWJlbCA9IHN5bWJvbCkpICsgCiAgICBnZW9tX2FibGluZShzbG9wZSA9IDEpICsKICAgIGdlb21fcG9pbnQoYWVzKGNvbG9yID0gZ3JvdXAsIHNpemUgPSBncm91cCkpICsgCiAgICBzY2FsZV9jb2xvcl9tYW51YWwoTlVMTCwgdmFsdWVzID0gc2NhdHRlci5jb2xvcnMpICsKICAgIHNjYWxlX3NpemVfbWFudWFsKHZhbHVlcyA9IHNjYXR0ZXIuc2l6ZSwgZ3VpZGUgPSAibm9uZSIpICsKICAgIGxhYnMoeCA9IGJxdW90ZShBW1BITzJdKSwgeSA9IGJxdW90ZShBW3BobzIqRGVsdGFdKSkgKwogICAgdGhlbWVfY293cGxvdCgpICsgcGFuZWxfYm9yZGVyKGNvbG9yID0gImdyYXkzMCIsIHNpemUgPSAxLjIpICsKICAgIHRoZW1lKGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMC44KSksCiAgICAgICAgICBsZWdlbmQucG9zaXRpb24gPSBjKDAuMDMsIDAuODMpLAogICAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChmYWNlID0gMiwgc2l6ZSA9IHJlbCgxLjIpKSwKICAgICAgICAgIGF4aXMubGluZSA9IGVsZW1lbnRfYmxhbmsoKSkKICByZXR1cm4ocCkKfQpgYGAKCmBgYHtyfQpwIDwtIG15X3NjYXR0ZXJfcGxvdF9maXgoKQpnZ3NhdmUoZmlsZW5hbWUgPSAiLi4vaW1nLzIwMjQwMzA4LWFsbC1jaGltZXJhLXNjYXR0ZXItY29sb3ItYnktUDJJRC5wbmciLAogICAgICAgcGxvdCA9IHAsIHdpZHRoID0gNC41LCBoZWlnaHQgPSA0LCBkcGkgPSAzMDApCmdncGxvdGx5KHAgKyBsYWJzKHggPSAiQTxzdWI+UEhPMjwvc3ViPiIsIHkgPSAiQTxzdWI+cGhvMjwvc3ViPiIpICsKICAgICAgICAgICB0aGVtZV9ncmF5KGJhc2Vfc2l6ZSA9IDE2KSArCiAgICAgICAgICAgdGhlbWUobGVnZW5kLnRleHQgPSBlbGVtZW50X21hcmtkb3duKCkpLCAKICAgICAgICAgdG9vbHRpcCA9IGMoImxhYmVsIiwgIngiLCAieSIpKQpgYGAKCnRoaXMgZnVuY3Rpb24gaXMgdGhlIHNhbWUgYXMgbXlfc2NhdHRlcl9wbG90X2ZpeCBleGNlcHQgdGhhdCBpdCBwbG90cyBhbGwgdGhlIGNoaW1lcmFzIHdpdGhvdXQgY29sb3JpbmcgdGhlbSBkaWZmZXJlbnRseS4gZm9yIGZpZ3VyZSA1CmBgYHtyfQpteV9zY2F0dGVyX3Bsb3RfYWxsIDwtIGZ1bmN0aW9uKCl7CiAgIyB0aGlzIGZ1bmN0aW9uIGlzIHRoZSBzYW1lIGFzIG15X3NjYXR0ZXJfcGxvdF9maXggZXhjZXB0IHRoYXQgaXQgcGxvdHMgYWxsIHRoZSBjaGltZXJhcwogICMgd2l0aG91dCBjb2xvcmluZyB0aGVtIGRpZmZlcmVudGx5LiBmb3IgZmlndXJlIDUKICBzMSA9IG15X2RhdGFfc2VsZWN0KHBhdHRlcm4gPSAiWFhYQ0NYIikKICBzMiA9IG15X2RhdGFfc2VsZWN0KHBhdHRlcm4gPSAiWFhYU1NYIikKICBzY2F0dGVyLmNvbG9ycyA9IGMoIlNjUGhvNCIgPSAiZm9yZXN0Z3JlZW4iLCAiQ2dQaG80IiA9ICJibHVlMyIsIAogICAgICAgICAgICAgICAgICAgICAiUDJJRDpDZyIgPSAiZ3JheTIwIiwgIlAySUQ6U2MiID0gImdyYXkyMCIsCiAgICAgICAgICAgICAgICAgICAgICJQMklEOm1peGVkIiA9ICJncmF5MjAiKQogIHNjYXR0ZXIuc2l6ZSA9IGMoIlNjUGhvNCIgPSAzLjUsICJDZ1BobzQiID0gMy41LAogICAgICAgICAgICAgICAgICAgIlAySUQ6Q2ciID0gMi41LCAiUDJJRDpTYyIgPSAyLjUsICJQMklEOm1peGVkIiA9IDIuNSkKICBwIDwtIHhpbWVyYSAlPiUgCiAgICBmaWx0ZXIoc2V0ICVpbiUgYygiTSIsICJTIikpICU+JSAKICAgIG11dGF0ZShBX1BITzIgPSBzaWduaWYoQV9QSE8yLCBkaWdpdHMgPSAyKSwKICAgICAgICAgICBBX3BobzIgPSBzaWduaWYoQV9waG8yLCBkaWdpdHMgPSAyKSwKICAgICAgICAgICBncm91cCA9IGNhc2Vfd2hlbigKICAgICAgICAgICAgIHN5bWJvbCA9PSAiQ0NDQ0MiIH4gIkNnUGhvNCIsCiAgICAgICAgICAgICBzeW1ib2wgPT0gIlNTU1NTIiB+ICJTY1BobzQiLAogICAgICAgICAgICAgcGxhc21pZCAlaW4lIHMxIH4gIlAySUQ6Q2ciLAogICAgICAgICAgICAgcGxhc21pZCAlaW4lIHMyIH4gIlAySUQ6U2MiLAogICAgICAgICAgICAgLmRlZmF1bHQgPSAiUDJJRDptaXhlZCIKICAgICAgICAgICApLAogICAgICAgICAgIGdyb3VwID0gZmN0X3JlbGV2ZWwoZ3JvdXAsIG5hbWVzKHNjYXR0ZXIuY29sb3JzKSkpICU+JSAKICAgIGdncGxvdChhZXMoeCA9IEFfUEhPMiwgeSA9IEFfcGhvMiwgbGFiZWwgPSBzeW1ib2wpKSArIAogICAgZ2VvbV9hYmxpbmUoc2xvcGUgPSAxKSArCiAgICBnZW9tX3BvaW50KGFlcyhjb2xvciA9IGdyb3VwLCBzaXplID0gZ3JvdXApKSArIAogICAgc2NhbGVfY29sb3JfbWFudWFsKE5VTEwsIHZhbHVlcyA9IHNjYXR0ZXIuY29sb3JzKSArCiAgICBzY2FsZV9zaXplX21hbnVhbCh2YWx1ZXMgPSBzY2F0dGVyLnNpemUsIGd1aWRlID0gIm5vbmUiKSArCiAgICBsYWJzKHggPSBicXVvdGUoQVtQSE8yXSksIHkgPSBicXVvdGUoQVtwaG8yXSkpICsKICAgIHRoZW1lX2Nvd3Bsb3QoKSArIHBhbmVsX2JvcmRlcihjb2xvciA9ICJncmF5MzAiLCBzaXplID0gMS4yKSArCiAgICB0aGVtZShsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDAuOCkpLAogICAgICAgICAgbGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiLAogICAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChmYWNlID0gMiwgc2l6ZSA9IHJlbCgxLjIpKSwKICAgICAgICAgIGF4aXMubGluZSA9IGVsZW1lbnRfYmxhbmsoKSkKCiAgcmV0dXJuKHApCn0KYGBgCgpQbG90IGFsbCBjaGltZXJhcywgZm9yIEZpZy4gNQpgYGB7cn0KcCA8LSBteV9zY2F0dGVyX3Bsb3RfYWxsKCkKZ2dzYXZlKGZpbGVuYW1lID0gIi4uL2ltZy8yMDI0MTEyMS1hbGwtY2hpbWVyYS1zY2F0dGVyLnBuZyIsCiAgICAgICBwbG90ID0gcCwgd2lkdGggPSA0LjUsIGhlaWdodCA9IDQsIGRwaSA9IDMwMCkKZ2dwbG90bHkocCArIGxhYnMoeCA9ICJBPHN1Yj5QSE8yPC9zdWI+IiwgeSA9ICJBPHN1Yj5waG8yPC9zdWI+IikgKwogICAgICAgICAgIHRoZW1lX2dyYXkoYmFzZV9zaXplID0gMTYpICsKICAgICAgICAgICB0aGVtZShsZWdlbmQudGV4dCA9IGVsZW1lbnRfbWFya2Rvd24oKSksIAogICAgICAgICB0b29sdGlwID0gYygibGFiZWwiLCAieCIsICJ5IikpCmBgYAojIyBTcG90bGlnaHQgaW5kaXZpZHVhbCBjaGltZXJhcwoKVGhlIGdvYWwgaGVyZSBpcyB0byBwbG90IGluZGl2aWR1YWwgY2hpbWVyYXMgaW4gb3JkZXIgdG8gdGVzdCBzcGVjaWZpYyBoeXBvdGhlc2VzIGFuZCBtYWtlIGNlcnRhaW4gcG9pbnRzLgoKMS4gIFdlIHNlcGFyYXRlbHkgdGVzdGVkIGFuZCBmb3VuZCB0aGF0IENnUGhvNCBEQkQgYmluZHMgdGhlIGNvbnNlbnN1cyBETkEgbW9yZSBzdHJvbmdseSB0aGFuIFNjUGhvNCBkb2VzLCBhbmQgaXQgYWxzbyBoYXMgdHdvIGFkZGl0aW9uYWwgYWN0aXZhdGlvbiBib29zdGVyIHJlZ2lvbnMsIHdoaWNoIGVuaGFuY2UgdGhlIGFjdGl2aXR5IG9mIHRoZSBtYWluIEFELiBXZSB0aGVyZWZvcmUgaHlwb3RoZXNpemUgdGhhdCBieSByZXBsYWNpbmcgdGhlIGNvcnJlc3BvbmRpbmcgcmVnaW9ucyBpbiBTY1BobzQgd2l0aCB0aGUgcGFydHMgZnJvbSBDZ1BobzQsIHdlIHdvdWxkIGNyZWF0ZSBhIGNoaW1lcmljIFRGIHRoYXQgaXMgbm90IG9yIGZhciBsZXNzIGRlcGVuZGVudCBvbiBQaG8yLgoyLiAgV2UgYWxzbyBleHBlY3QgdGhhdCB0aG9zZSByZWdpb25zIGFkZGl0aXZlbHkgY29udHJpYnV0ZSB0byB0aGUgcmVkdWNlZCBQaG8yLWRlcGVuZGVuY2UsIHNob3duIGFzIGluY3JlYXNlZCBURiBhY3Rpdml0eSBvZiB0aGUgY2hpbWVyYSBpbiB0aGUgKnBobzLiiIYqIGJhY2tncm91bmQuCgpEZXNpZ24gcGxvdAoKYGBge3J9Cm15X3Bsb3Rfc3Vic2V0X3hpbWVyYSA8LSBmdW5jdGlvbihzeW1ib2xzKXsKICAjIHRoaXMgZnVuY3Rpb24gcGxvdHMgYSBzdWJzZXQgb2YgdGhlIGNoaW1lcmFzIGFzIGhvcml6b250YWwgYmFyIHBsb3RzCiAgIyBzaG93aW5nIHRoZSBSZWwuIEFfUEhPMiBhbmQgJUFfcGhvMuKIhiB2YWx1ZXMKICAjIGl0IHRha2VzIGFzIGlucHV0IGEgdmVjdG9yIGNvbnRhaW5pbmcgdGhlIHN5bWJvbHMgZm9yIHRoZSBjaGltZXJhcyBmb3IgCiAgIyBwbG90dGluZy4gdGhlIG9yZGVyIGluIHRoZSB2ZWN0b3IgZGV0ZXJtaW5lcyB0aGUgcGxvdCBvcmRlcgogICMgdGhlIGVuZG9nZW5vdXMgU2NQaG80IGFuZCBDZ1BobzQgYXJlIGltcGxpZWQKICBtaXNzaW5nIDwtIHNldGRpZmYoc3ltYm9scywgeGltZXJhJHN5bWJvbCkKICBpZihsZW5ndGgobWlzc2luZykgIT0gMCkKICAgIHN0b3AocGFzdGUobWlzc2luZywgImFyZSBub3QgZm91bmQiLCBzZXAgPSAiICIpKQogIAogIHRtcCA8LSBmaWx0ZXIoeGltZXJhLCBzeW1ib2wgJWluJSBjKCJTU1NTUyIsICJDQ0NDQyIsIHN5bWJvbHMpKSAlPiUgCiAgICBtdXRhdGUoCiAgICAgIHJTRV9QSE8yID0gc2VfUEhPMiAvIEFfUEhPMltzeW1ib2wgPT0gIlNTU1NTIl0sCiAgICAgIHJTRV9waG8yID0gc2VfcGhvMiAvIEFfcGhvMltzeW1ib2wgPT0gIlNTU1NTIl0KICAgICkgJT4lIAogICAgcGl2b3RfbG9uZ2VyKGNvbHMgPSBjKHJBX1BITzIsIHJBX3BobzIsIHJTRV9QSE8yLCByU0VfcGhvMiksIAogICAgICAgICAgICAgICAgICNwaXZvdF9sb25nZXIoY29scyA9IGMoQV9QSE8yLCBBX3BobzIsIHNlX1BITzIsIHNlX3BobzIpLCAKICAgICAgICAgICAgICAgICBuYW1lc190byA9IGMoIi52YWx1ZSIsICJwYXJhbWV0ZXIiKSwgbmFtZXNfc2VwID0gIl8iLAogICAgICAgICAgICAgICAgIHZhbHVlc190byA9ICJ2YWx1ZSIpICU+JSAKICAgIG11dGF0ZShwYXJhbWV0ZXIgPSBmY3RfcmVsZXZlbChwYXJhbWV0ZXIsICJQSE8yIiksCiAgICAgICAgICAgc3ltYm9sID0gZmFjdG9yKHN5bWJvbCwgbGV2ZWxzID0gCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5pcXVlKGMoIlNTU1NTIiwgIkNDQ0NDIiwgc3ltYm9scykpKSkgJT4lIAogICAgc2VsZWN0KC1jKEFfUEhPMjpib29zdCkpCiAgCiAgIyBsYWJlbGxlcgogIHBhci5leHBsYWluIDwtIGMoCiAgICBQSE8yID0gIlJlbC4gQTxzdWI+UEhPMjwvc3ViPiIsCiAgICAjYm9vc3QgPSAiQm9vc3QiLAogICAgcGhvMiA9ICJSZWwuIEE8c3ViPnBobzLiiIY8L3N1Yj4iCiAgKQogIAogIHAgPC0gZ2dwbG90KHRtcCwgYWVzKHkgPSBzeW1ib2wsIHggPSByQSkpICsKICAgIGdlb21fY29sKHdpZHRoID0gMC41LCBjb2xvciA9ICJibGFjayIsIGZpbGwgPSAiZ3JheTgwIikgKwogICAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMSwgbGluZXR5cGUgPSAyLCBjb2xvciA9ICJncmF5MzAiKSArCiAgICBnZW9tX2Vycm9yYmFyKGFlcyh4bWluID0gckEgLSByU0UsIHhtYXggPSByQSArIHJTRSksIHdpZHRoID0gMC4yKSArCiAgICBmYWNldF93cmFwKH5wYXJhbWV0ZXIsIHNjYWxlcyA9ICJmcmVlX3giLCMgc3dpdGNoID0gIngiLAogICAgICAgICAgICAgIGxhYmVsbGVyID0gbGFiZWxsZXIocGFyYW1ldGVyID0gcGFyLmV4cGxhaW4pKSArCiAgICBzY2FsZV95X2Rpc2NyZXRlKGxpbWl0cyA9IHJldikgKyAKICAgIHNjYWxlX3hfY29udGludW91cyhleHBhbmQgPSBleHBhbnNpb24obXVsdCA9IGMoMC4wMiwgMC4wNSkpKSArCiAgICB0aGVtZV9jb3dwbG90KCkgKyBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIikgKwogICAgYmFja2dyb3VuZF9ncmlkKG1ham9yID0gInkiLCBtaW5vciA9ICJub25lIikgKwogICAgdGhlbWUoYXhpcy50ZXh0LnkgPSBlbGVtZW50X3RleHQoZmFtaWx5ID0gImNvdXJpZXIiKSwKICAgICAgICAgIGF4aXMudGl0bGUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICBzdHJpcC5wbGFjZW1lbnQgPSAib3V0c2lkZSIsCiAgICAgICAgICBzdHJpcC5iYWNrZ3JvdW5kID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgICAgc3RyaXAudGV4dCA9IGVsZW1lbnRfbWFya2Rvd24oKSkKICByZXR1cm4ocCkKfQpgYGAKCjxmb250IGNvbG9yPSJyZWQiPlVwZGF0ZSAyMDI0LTExLTIyPC9mb250PgoKQWx0ZXJuYXRpdmUgZGVzaWduLCB3aXRoIGluZGl2aWR1YWwgcG9pbnRzIGFuZCBub3QgcmVsYXRpdmUgdG8gU2NQaG80LiBBbHNvLCBpbmRpdmlkdWFsIGRhdGFwb2ludHMgd2VyZSBwbG90dGVkIGZvciBzYW1wbGVzIHdpdGggXDwxMCByZXBsaWNhdGVzCgpgYGB7cn0KbXlfcGxvdF9zdWJzZXRfeGltZXJhX2FsdCA8LSBmdW5jdGlvbihzeW1ib2xzKXsKICAjIHRoaXMgZnVuY3Rpb24gcGxvdHMgYSBzdWJzZXQgb2YgdGhlIGNoaW1lcmFzIGFzIGhvcml6b250YWwgYmFyIHBsb3RzCiAgIyBzaG93aW5nIHRoZSBSZWwuIEFfUEhPMiBhbmQgJUFfcGhvMuKIhiB2YWx1ZXMKICAjIGl0IHRha2VzIGFzIGlucHV0IGEgdmVjdG9yIGNvbnRhaW5pbmcgdGhlIHN5bWJvbHMgZm9yIHRoZSBjaGltZXJhcyBmb3IgCiAgIyBwbG90dGluZy4gdGhlIG9yZGVyIGluIHRoZSB2ZWN0b3IgZGV0ZXJtaW5lcyB0aGUgcGxvdCBvcmRlcgogICMgdGhlIGVuZG9nZW5vdXMgU2NQaG80IGFuZCBDZ1BobzQgYXJlIGltcGxpZWQKICBtaXNzaW5nIDwtIHNldGRpZmYoc3ltYm9scywgeGltZXJhJHN5bWJvbCkKICBpZihsZW5ndGgobWlzc2luZykgIT0gMCkKICAgIHN0b3AocGFzdGUobWlzc2luZywgImFyZSBub3QgZm91bmQiLCBzZXAgPSAiICIpKQogIAogIHRtcCA8LSBmaWx0ZXIoZGF0X3NlcCwgc3ltYm9sICVpbiUgYygiU1NTU1MiLCAiQ0NDQ0MiLCBzeW1ib2xzKSkgJT4lIAogICAgbXV0YXRlKGhvc3QgPSBmY3RfcmVsZXZlbChob3N0LCAiUEhPMiIpLAogICAgICAgICAgIHN5bWJvbCA9IGZhY3RvcihzeW1ib2wsIGxldmVscyA9IAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVuaXF1ZShjKCJTU1NTUyIsICJDQ0NDQyIsIHN5bWJvbHMpKSkpCgogIHRtcCAlPiUgY291bnQoc3ltYm9sLCBob3N0KSAlPiUgcHJpbnQoKQogICMgbGFiZWxsZXIKICBwYXIuZXhwbGFpbiA8LSBjKAogICAgUEhPMiA9ICJBPHN1Yj5QSE8yPC9zdWI+IiwKICAgICNib29zdCA9ICJCb29zdCIsCiAgICBwaG8yID0gIkE8c3ViPnBobzLiiIY8L3N1Yj4iCiAgKQogIAogIHAgPC0gZ2dwbG90KHRtcCwgYWVzKHkgPSBzeW1ib2wsIHggPSBBKSkgKwogICAgZ2VvbV9iYXIoc3RhdCA9ICJzdW1tYXJ5IiwgZnVuID0gIm1lYW4iLCAKICAgICAgICAgICAgIHdpZHRoID0gMC41LCBjb2xvciA9ICJibGFjayIsIGZpbGwgPSAiZ3JheTgwIikgKwogICAgc3RhdF9zdW1tYXJ5KGZ1bi5kYXRhID0gIm1lYW5fY2xfYm9vdCIsIGdlb20gPSAibGluZXJhbmdlIiwKICAgICAgICAgICAgICAgICBjb2xvciA9ICJzdGVlbGJsdWU0IikgKwogICAgZ2VvbV9wb2ludChkYXRhID0gZmlsdGVyKHRtcCwgIXN5bWJvbCAlaW4lIGMoIkNDQ0NDIiwgIlNTU1NTIikpLCAKICAgICAgICAgICAgICAgc2l6ZSA9IDAuNiwgc2hhcGUgPSAzLCBjb2xvciA9ICJncmF5MzAiKSArCiAgICAjZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMSwgbGluZXR5cGUgPSAyLCBjb2xvciA9ICJncmF5MzAiKSArCiAgICAjZ2VvbV9lcnJvcmJhcihhZXMoeG1pbiA9IHJBIC0gclNFLCB4bWF4ID0gckEgKyByU0UpLCB3aWR0aCA9IDAuMikgKwogICAgZmFjZXRfd3JhcCh+aG9zdCwgc2NhbGVzID0gImZyZWVfeCIsIyBzd2l0Y2ggPSAieCIsCiAgICAgICAgICAgICAgbGFiZWxsZXIgPSBsYWJlbGxlcihob3N0ID0gcGFyLmV4cGxhaW4pKSArCiAgICBzY2FsZV95X2Rpc2NyZXRlKGxpbWl0cyA9IHJldikgKyAKICAgIHNjYWxlX3hfY29udGludW91cyhleHBhbmQgPSBleHBhbnNpb24obXVsdCA9IGMoMC4wMiwgMC4wNSkpKSArCiAgICB0aGVtZV9jb3dwbG90KCkgKyBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIikgKwogICAgYmFja2dyb3VuZF9ncmlkKG1ham9yID0gInkiLCBtaW5vciA9ICJub25lIikgKwogICAgdGhlbWUoYXhpcy50ZXh0LnkgPSBlbGVtZW50X3RleHQoZmFtaWx5ID0gImNvdXJpZXIiKSwKICAgICAgICAgIGF4aXMudGl0bGUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICBzdHJpcC5wbGFjZW1lbnQgPSAib3V0c2lkZSIsCiAgICAgICAgICBzdHJpcC5iYWNrZ3JvdW5kID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgICAgc3RyaXAudGV4dCA9IGVsZW1lbnRfbWFya2Rvd24oKSkKICByZXR1cm4ocCkKfQpgYGAKCiMjIyBNaW5pbWFsIENnUGhvNCBwYXJ0cyBmb3IgQV9waG8yCgpUaGUgY2hpbWVyYSB3aXRoIHRoZSBsZWFzdCBhbW91bnQgb2YgQ2dQaG80IGFuZCB5ZXQgaGF2ZSBhcHByZWNpYWJsZSBhY3Rpdml0eSBpbiB0aGUgYWJzZW5jZSBvZiBQaG8yIGlzIFRoZXNlIGluY2x1ZGUgU1NTU1MsIENDQ0NDLCBTU1NTQywgQ1NTU1MsIFNTQ1NTLCBDU0NTUywgQ1NTU0MsIENTQ1NDCgpgYGB7cn0Kc2VsZWN0ZWQgPC0gYXMuY2hhcmFjdGVyKAogIGV4cHJlc3Npb24oQ1NTU1MsIFNDU1NTLCBDQ1NTUywgU1NDU1MsIENTQ1NTLCBTQ0NTUywgQ0NDU1MsIFNTU1NDLCBTU0NTQywgQ1NDU0MsIENTU2NDKSkKI3NlbGVjdGVkIDwtIGZpbHRlcihtZXRhLCBzeW1ib2wgJWluJSBzZWxlY3RlZCkgJT4lIHB1bGwocGxhc21pZCkKbXlfcGxvdF9zdWJzZXRfeGltZXJhKHNlbGVjdGVkKQpnZ3NhdmUoIi4uL2ltZy8yMDI0MDMwOC1zZWxlY3RlZC1jaGltZXJhLXJlbC1hY3Rpdml0eS5wbmciLCB3aWR0aCA9IDQsIGhlaWdodCA9IDQpCiMgcGxvdCB3aXRoIGFic29sdXRlIEEgbm90IHJlbGF0aXZlLCBhbmQgcGxvdCBpbmRpdmlkdWFsIGRhdGEgcG9pbnRzCm15X3Bsb3Rfc3Vic2V0X3hpbWVyYV9hbHQoc2VsZWN0ZWQpCmdnc2F2ZSgiLi4vaW1nLzIwMjQxMTIyLXNlbGVjdGVkLWNoaW1lcmEtcmVsLWFjdGl2aXR5LnBuZyIsIHdpZHRoID0gNCwgaGVpZ2h0ID0gNCkKYGBgCgojIyMgUmVnaW9uIDEtMyBtYWluIGVmZmVjdHMgYW5kIGludGVyYWN0aW9ucwoKYGBge3J9CiMgc2VsZWN0IHRoZSBjaGltZXJhcwpzZWxlY3RlZCA8LSBhcy5jaGFyYWN0ZXIoCiAgZXhwcmVzc2lvbihTU1NTUywgQ1NTU1MsIFNDU1NTLCBTU0NTUywgQ0NTU1MsIENTQ1NTLCBTQ0NTUywgQ0NDU1MpCikKCiMgZXh0cmFjdCB0aGUgZGF0YQp0bXAgPC0geGltZXJhICU+JSAKICBmaWx0ZXIoc3ltYm9sICVpbiUgc2VsZWN0ZWQsIHNldCA9PSAiTSIpICU+JSAKICBzZWxlY3QocGxhc21pZCwgc3ltYm9sLCBncm91cCkgJT4lIAogIGlubmVyX2pvaW4oZGF0LCBieSA9ICJwbGFzbWlkIikgJT4lIAogIG11dGF0ZSggYFIvR2AgPSBZTDIuSCAvIEJMMS5IICkgJT4lCiAgZmlsdGVyKGZsYWcgPT0gInBhc3MiKSAlPiUgCiAgc2VsZWN0KC1uUkZQLCAtbkdGUCwgLXdlbGwsIC1mbGFnKQoKIyBwcmVwYXJlIHRoZSBmYWN0b3IgbGV2ZWxzCnNwbGl0IDwtIGMoMSwxLDEsMik7IG5hbWVzKHNwbGl0KSA8LSBjKCJSMSIsICJBRCIsICJOTFMiKQoKdG1wIDwtIHRtcCAlPiUgCiAgc2VwYXJhdGVfd2lkZXJfcG9zaXRpb24oc3ltYm9sLCBzcGxpdCkgJT4lIAogIG11dGF0ZShhY3Jvc3MoUjE6TkxTLCB+ZmFjdG9yKC54LCBsZXZlbHMgPSBjKCJTIiwgIkMiKSkpKQoKIyB0ZXN0IEFfUEhPMgpwcmludCgiVGVzdGluZyBBX1BITzIiKQpsbS5yZXMgPC0gdG1wICU+JSAKICBmaWx0ZXIoaG9zdCA9PSAiUEhPMiIpICU+JSAKICBsbShgUi9HYCB+IChSMSpBRCpOTFMpLCBkYXRhID0gLikgJT4lIAogIHN1bW1hcnkoKQojIGFkZGluZyBhZGp1c3RlZCBQLXZhbHVlCmxtLnJlcyRjb2VmZmljaWVudHMgPC0gY2JpbmQoCiAgY29lZihsbS5yZXMpLAogICJQLmFkaiIgPSBwLmFkanVzdChjb2VmKGxtLnJlcylbLCdQcig+fHR8KSddLCBtZXRob2QgPSAiaG9sbSIpCikKcHJpbnQobG0ucmVzKQojIHN0b3JlIHRoZSB0ZXN0IHJlc3VsdHMgZm9yIHBsb3R0aW5nCnJlcy5QSE8yIDwtIGNvZWYobG0ucmVzKVstMSxdICU+JSBhc190aWJibGUocm93bmFtZXMgPSAiY29tcG9uZW50IikKCiMgdGVzdCBBX3BobzIKcHJpbnQoIlRlc3RpbmcgQV9waG8y4oiGIikKbG0ucmVzIDwtIHRtcCAlPiUgCiAgZmlsdGVyKGhvc3QgPT0gInBobzIiKSAlPiUgCiAgbG0oYFIvR2AgfiAoUjEqQUQqTkxTKSwgZGF0YSA9IC4pICU+JSAKICBzdW1tYXJ5KCkKIyBhZGRpbmcgYWRqdXN0ZWQgUC12YWx1ZQpsbS5yZXMkY29lZmZpY2llbnRzIDwtIGNiaW5kKAogIGNvZWYobG0ucmVzKSwKICAiUC5hZGoiID0gcC5hZGp1c3QoY29lZihsbS5yZXMpWywnUHIoPnx0fCknXSwgbWV0aG9kID0gImhvbG0iKQopCnByaW50KGxtLnJlcykKIyBzdG9yZSB0aGUgdGVzdCByZXN1bHRzIGZvciBwbG90dGluZwpyZXMucGhvMiA8LSBjb2VmKGxtLnJlcylbLTEsXSAlPiUgYXNfdGliYmxlKHJvd25hbWVzID0gImNvbXBvbmVudCIpCgojIGNvbWJpbmUgdGhlIHJlc3VsdHMKdGVzdC5yZXMgPC0gYmluZF9yb3dzKAogICJBX1BITzIiID0gcmVzLlBITzIsICJBX3BobzIiID0gcmVzLnBobzIsIC5pZCA9ICJwYXJhbWV0ZXIiCikKYGBgCgpQbG90IHRoZSByZXN1bHQKCmBgYHtyfQpwYXIuZXhwbGFpbiA8LSBjKAogIEFfUEhPMiA9ICJBPHN1Yj5QSE8yPC9zdWI+IiwKICBBX3BobzIgPSAiQTxzdWI+cGhvMuKIhjwvc3ViPiIKKQoKcCA8LSB0ZXN0LnJlcyAlPiUgCiAgcmVuYW1lKGVzdGltYXRlID0gRXN0aW1hdGUsIHNlID0gYFN0ZC4gRXJyb3JgKSAlPiUgCiAgbXV0YXRlKAogICAgY29tcG9uZW50ID0gZ3N1YigiQyIsICIiLCBjb21wb25lbnQpICU+JSBmY3RfaW5vcmRlcigpLAogICAgcGFyYW1ldGVyID0gZmFjdG9yKHBhcmFtZXRlciwgbGV2ZWxzID0gYygiQV9QSE8yIiwgIkFfcGhvMiIpKSwKICAgIHNpZyA9IFAuYWRqIDwgMC4wNQogICkgJT4lIAogIGdncGxvdChhZXMoeCA9IGNvbXBvbmVudCwgeSA9IGVzdGltYXRlKSkgKwogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDAsIGxpbmV0eXBlID0gMSwgY29sb3IgPSAiZ3JheTUwIikgKwogIGdlb21fY29sKGFlcyhmaWxsID0gUC5hZGogPCAwLjA1KSwgd2lkdGggPSAwLjUsIGNvbG9yID0gImJsYWNrIikgKwogIGdlb21fcG9pbnRyYW5nZShhZXMoeW1pbiA9IGVzdGltYXRlLXNlLCB5bWF4ID0gZXN0aW1hdGUrc2UpLCBzaXplID0gMC4yKSArCiAgZmFjZXRfd3JhcCh+cGFyYW1ldGVyLCBzY2FsZXMgPSAiZnJlZV95IiwgbnJvdyA9IDIsCiAgICAgICAgICAgICBsYWJlbGxlciA9IGxhYmVsbGVyKHBhcmFtZXRlciA9IHBhci5leHBsYWluKSkgKwogIHNjYWxlX3hfZGlzY3JldGUoKSArIAogIHNjYWxlX3lfY29udGludW91cygpICsKICBzY2FsZV9maWxsX21hbnVhbChOVUxMLCAKICAgICAgICAgICAgICAgICAgICB2YWx1ZXMgPSBjKCJncmF5OTAiLCAiZ3JheTUwIikpICsKICB0aGVtZV9jb3dwbG90KCkgKyBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIikgKwogIGJhY2tncm91bmRfZ3JpZChtYWpvciA9ICJ5IiwgbWlub3IgPSAieSIpICsKICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCBoanVzdCA9IDEsIHNpemUgPSByZWwoMSkpLAogICAgICAgIGF4aXMudGl0bGUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgYXhpcy5saW5lID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgIGxlZ2VuZC5wb3NpdGlvbiA9ICJib3R0b20iLAogICAgICAgIHN0cmlwLnBsYWNlbWVudCA9ICJvdXRzaWRlIiwKICAgICAgICBzdHJpcC5iYWNrZ3JvdW5kID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgIHN0cmlwLnRleHQgPSBlbGVtZW50X21hcmtkb3duKHNpemUgPSByZWwoMSkpKQpwCmdnc2F2ZSgiLi4vaW1nLzIwMjQwNDA0LXJlZ2lvbjEtMy1lcGlzdGFzaXMtcGxvdC5wbmciLCB3aWR0aCA9IDMuNSwgaGVpZ2h0ID0gNC41KQpgYGAKCiMjIyBSZWdpb24gNCBzcGxpdHMKCldlIGhhdmUgc28gZmFyIGZvY3VzZWQgb24gdGhlIG1haW4gc2V0IHdpdGggdGhlIDUgcmVnaW9uIGRlc2lnbi4gSW4gdGhlIHNjYXR0ZXIgcGxvdCBiZWxvdywgd2Ugc2VlIHRoYXQgdGhlcmUgaXMgYSBzdWJzZXQgb2YgY2hpbWVyYXMgaW4gYmV0d2VlbiB0aGUgUDJJRDpTYyBhbmQgUDJJRDpDZyBvbmVzLiBUaGV5IGFyZSBpbnRlcmVzdGluZyBpbiB0aGF0IHRoZWlyIEFfcGhvMuKIhi9BX1BITzIgcmF0aW9zIGFyZSBpbnRlcm1lZGlhdGUuCgohW3NjYXR0ZXJdKC4uL2ltZy8yMDIzMTIyMC1hbGwtY2hpbWVyYS1zY2F0dGVyLWNvbG9yLWJ5LVAySUQucG5nKQoKYGBge3J9CnNlbGVjdGVkIDwtIGFzLmNoYXJhY3RlcigKICBleHByZXNzaW9uKENDQ1NDLCBDQ0Njc0MsIENDQ3NjQywKICAgIENTU0NDLCBDU1NTQywgQ1NTY3NDLCBDU1NzY0MKICAgICAgICAgICAgICkpCiNzZWxlY3RlZCA8LSBmaWx0ZXIobWV0YSwgc3ltYm9sICVpbiUgc2VsZWN0ZWQpICU+JSBwdWxsKHBsYXNtaWQpCm15X3Bsb3Rfc3Vic2V0X3hpbWVyYV9hbHQoc2VsZWN0ZWQpICsgc2NhbGVfeF9jb250aW51b3VzKGV4cGFuZCA9IGV4cGFuc2lvbihtdWx0ID0gYygwLjAyLCAwLjE1KSkpCmdnc2F2ZSgiLi4vaW1nLzIwMjQxMTIyLVAySUQtc3BsaXQtcmVsLWFjdGl2aXR5LnBuZyIsIHdpZHRoID0gMy41LCBoZWlnaHQgPSAzLjUpCmBgYAoKU3RhdGlzdGljYWwgdGVzdHMgZm9yIGdyb3VwIDEKCmBgYHtyfQojIHNlbGVjdCB0aGUgY2hpbWVyYXMKc2VsZWN0ZWQgPC0gYXMuY2hhcmFjdGVyKAogIGV4cHJlc3Npb24oQ0NDQ0MsIENDQ1NDLCBDQ0Njc0MsIENDQ3NjQykKKQoKIyBleHRyYWN0IHRoZSBkYXRhCnRtcCA8LSB4aW1lcmEgJT4lIAogIGZpbHRlcihzeW1ib2wgJWluJSBzZWxlY3RlZCkgJT4lIAogIHNlbGVjdChwbGFzbWlkLCBzeW1ib2wsIGdyb3VwKSAlPiUgCiAgaW5uZXJfam9pbihkYXQsIGJ5ID0gInBsYXNtaWQiKSAlPiUgCiAgbXV0YXRlKCBgUi9HYCA9IFlMMi5IIC8gQkwxLkggKSAlPiUKICBmaWx0ZXIoZmxhZyA9PSAicGFzcyIpICU+JSAKICBzZWxlY3QoLW5SRlAsIC1uR0ZQLCAtd2VsbCwgLWZsYWcpICU+JSAKICBtdXRhdGUoc3ltYm9sID0gZmFjdG9yKHN5bWJvbCwgbGV2ZWxzID0gISFzZWxlY3RlZCkpCgojIHRlc3QgQV9QSE8yCnByaW50KCJUZXN0aW5nIEFfUEhPMiIpCmxtLnJlcyA8LSB0bXAgJT4lIAogIGZpbHRlcihob3N0ID09ICJQSE8yIikgJT4lIAogIGxtKGBSL0dgIH4gc3ltYm9sLCBkYXRhID0gLikgJT4lIAogIHN1bW1hcnkoKQojIGFkZGluZyBhZGp1c3RlZCBQLXZhbHVlCmxtLnJlcyRjb2VmZmljaWVudHMgPC0gY2JpbmQoCiAgY29lZihsbS5yZXMpLAogICJQLmFkaiIgPSBwLmFkanVzdChjb2VmKGxtLnJlcylbLCdQcig+fHR8KSddLCBtZXRob2QgPSAiaG9sbSIpCikKcHJpbnQobG0ucmVzKQojIHN0b3JlIHRoZSB0ZXN0IHJlc3VsdHMgZm9yIHBsb3R0aW5nCiNyZXMuUEhPMiA8LSBjb2VmKGxtLnJlcylbLTEsXSAlPiUgYXNfdGliYmxlKHJvd25hbWVzID0gImNvbXBvbmVudCIpCgojIHRlc3QgQV9waG8yCnByaW50KCJUZXN0aW5nIEFfcGhvMuKIhiIpCmxtLnJlcyA8LSB0bXAgJT4lIAogIGZpbHRlcihob3N0ID09ICJwaG8yIikgJT4lIAogIGxtKGBSL0dgIH4gc3ltYm9sLCBkYXRhID0gLikgJT4lIAogIHN1bW1hcnkoKQojIGFkZGluZyBhZGp1c3RlZCBQLXZhbHVlCmxtLnJlcyRjb2VmZmljaWVudHMgPC0gY2JpbmQoCiAgY29lZihsbS5yZXMpLAogICJQLmFkaiIgPSBwLmFkanVzdChjb2VmKGxtLnJlcylbLCdQcig+fHR8KSddLCBtZXRob2QgPSAiaG9sbSIpCikKcHJpbnQobG0ucmVzKQpgYGAKClN0YXRpc3RpY2FsIHRlc3RzIGZvciBncm91cCAyCgpgYGB7cn0KIyBzZWxlY3QgdGhlIGNoaW1lcmFzCnNlbGVjdGVkIDwtIGFzLmNoYXJhY3RlcigKICBleHByZXNzaW9uKCBDU1NDQywgQ1NTU0MsIENTU2NzQywgQ1NTc2NDICkKKQoKIyBleHRyYWN0IHRoZSBkYXRhCnRtcCA8LSB4aW1lcmEgJT4lIAogIGZpbHRlcihzeW1ib2wgJWluJSBzZWxlY3RlZCkgJT4lIAogIHNlbGVjdChwbGFzbWlkLCBzeW1ib2wsIGdyb3VwKSAlPiUgCiAgaW5uZXJfam9pbihkYXQsIGJ5ID0gInBsYXNtaWQiKSAlPiUgCiAgbXV0YXRlKCBgUi9HYCA9IFlMMi5IIC8gQkwxLkggKSAlPiUKICBmaWx0ZXIoZmxhZyA9PSAicGFzcyIpICU+JSAKICBzZWxlY3QoLW5SRlAsIC1uR0ZQLCAtd2VsbCwgLWZsYWcpICU+JSAKICBtdXRhdGUoc3ltYm9sID0gZmFjdG9yKHN5bWJvbCwgbGV2ZWxzID0gISFzZWxlY3RlZCkpCgojIHRlc3QgQV9QSE8yCnByaW50KCJUZXN0aW5nIEFfUEhPMiIpCmxtLnJlcyA8LSB0bXAgJT4lIAogIGZpbHRlcihob3N0ID09ICJQSE8yIikgJT4lIAogIGxtKGBSL0dgIH4gc3ltYm9sLCBkYXRhID0gLikgJT4lIAogIHN1bW1hcnkoKQojIGFkZGluZyBhZGp1c3RlZCBQLXZhbHVlCmxtLnJlcyRjb2VmZmljaWVudHMgPC0gY2JpbmQoCiAgY29lZihsbS5yZXMpLAogICJQLmFkaiIgPSBwLmFkanVzdChjb2VmKGxtLnJlcylbLCdQcig+fHR8KSddLCBtZXRob2QgPSAiaG9sbSIpCikKcHJpbnQobG0ucmVzKQojIHN0b3JlIHRoZSB0ZXN0IHJlc3VsdHMgZm9yIHBsb3R0aW5nCiNyZXMuUEhPMiA8LSBjb2VmKGxtLnJlcylbLTEsXSAlPiUgYXNfdGliYmxlKHJvd25hbWVzID0gImNvbXBvbmVudCIpCgojIHRlc3QgQV9waG8yCnByaW50KCJUZXN0aW5nIEFfcGhvMuKIhiIpCmxtLnJlcyA8LSB0bXAgJT4lIAogIGZpbHRlcihob3N0ID09ICJwaG8yIikgJT4lIAogIGxtKGBSL0dgIH4gc3ltYm9sLCBkYXRhID0gLikgJT4lIAogIHN1bW1hcnkoKQojIGFkZGluZyBhZGp1c3RlZCBQLXZhbHVlCmxtLnJlcyRjb2VmZmljaWVudHMgPC0gY2JpbmQoCiAgY29lZihsbS5yZXMpLAogICJQLmFkaiIgPSBwLmFkanVzdChjb2VmKGxtLnJlcylbLCdQcig+fHR8KSddLCBtZXRob2QgPSAiaG9sbSIpCikKcHJpbnQobG0ucmVzKQpgYGAKCiMjIFJlZ2lvbiBtYWluIGVmZmVjdAoKYGBge3J9CnNwbGl0IDwtIGMoMSwxLDEsMSwxKTsgbmFtZXMoc3BsaXQpIDwtIHBhc3RlMCgiUCIsIDE6NSkKdG1wIDwtIHhpbWVyYSAlPiUgCiAgZmlsdGVyKHNldCA9PSAiTSIsIGdyb3VwICE9ICJuLmYuIikgJT4lIAogIHNlcGFyYXRlX3dpZGVyX3Bvc2l0aW9uKHN5bWJvbCwgc3BsaXQpICU+JSAKICBtdXRhdGUoYWNyb3NzKFAxOlA1LCB+ZmFjdG9yKC54LCBsZXZlbHMgPSBjKCJTIiwgIkMiKSkpKQpsbSA8LSBsbShBX3BobzIgfiAoUDErUDIrUDMrUDQrUDUpLCBkYXRhID0gdG1wKQpzdW1tYXJ5KGxtKQpgYGAKClRoZSBtYWluIGVmZmVjdHMgd2VyZSBjYWxjdWxhdGVkIGJ5IGF2ZXJhZ2luZyBvdmVyIGFsbCBjaGltZXJhcyB3aXRoIENnUGhvNCByZWdpb24gYXQgdGhlIHJlc3BlY3RpdmUgcG9zaXRpb24uIEknZCBsaWtlIHRvIGJyZWFrIHRoZW0gZG93biBieSBiYWNrZ3JvdW5kcy4gRm9yIGV4YW1wbGUsIGZvciByZWdpb24gMywgSSdkIGxpa2UgdG8gc2VlIHRoZSBwYWlyd2lzZSBjb21wYXJpc29ucyBiZXR3ZWVuIENDQ1NTIGFuZCBDQ1NTUywgd2hlcmUgb25seSByZWdpb24gMyBkaWZmZXJzLiBUaGUgc3RlcHMgYXJlCgoxLiAgc2VsZWN0IHRoZSByZWdpb24gdG8gYmUgY29tcGFyZWQuIHNwbGl0IHRoZSBzeW1ib2wgaW50byB0d28gcGFydHMgLSB0aGUgZ2Vub3R5cGUgb2YgdGhlIGZvY2FsIHJlZ2lvbiBhbmQgdGhlIHJlc3QKMi4gIGdyb3VwIGJ5IHRoZSBzZWNvbmQgcGFydCAocmVzdCkgYW5kIGNhbGN1bGF0ZSB0aGUgZGlmZmVyZW50aWFsCgpgYGB7cn0KbXlfY2FsY19yZWdpb25fZWZmZWN0IDwtIGZ1bmN0aW9uKHJlZ2lvbiwgdmFyaWFibGUpewogICMgdGhpcyBmdW5jdGlvbiB0YWtlcyB0aGUgbmFtZSBvZiBhIHZhcmlhYmxlIG9mIGludGVyZXN0CiAgIyB4IHNwZWNpZmllcyB0aGUgZm9yZWdyb3VuZCByZWdpb24sIHdoaWNoIHdpbGwgYmUgZXhhbWluZWQgZm9yIGl0cyBlZmZlY3Qgb24KICAjIHRoZSB2YXJpYWJsZSBvZiBpbnRlcmVzdC4KICAjIGl0IHRoZW4gdHJhbnNmb3JtcyB0aGUgeGltZXJhIGRhdGEgZnJhbWUgdG8gcHJlc2VydmUgb25seSB0aGUgdmFyaWFibGUgb2YKICAjIGludGVyZXN0LCBwaXZvdHMgaXQgd2lkZXIgYWZ0ZXIgZ3JvdXBpbmcgYnkgdGhlIGJhY2tncm91bmQgY29tcG9zaXRpb24uCiAgCiAgIyBwcmVwYXJlIHRoZSBkYXRhIGJ5IG11dGF0aW5nIHRoZSBzeW1ib2wgY29sdW1uIGludG8gZmcgYW5kIGJnCiAgdmFsaWQudmFyIDwtIGMoIkFfUEhPMiIsICJBX3BobzIiLCAickFfUEhPMiIsICJyQV9waG8yIiwgImJvb3N0IikKICBpZighdmFyaWFibGUgJWluJSB2YWxpZC52YXIpCiAgICBzdG9wKHBhc3RlMCgiUGxlYXNlIHNwZWNpZnkgb25lIG9mIHRoZSB2YWxpZCB2YXJpYWJsZSBuYW1lczoiLCAKICAgICAgICAgICAgICAgIHBhc3RlKHZhbGlkLnZhciwgY29sbGFwc2UgPSAiLCAiKSkpCiAgdG1wIDwtIHhpbWVyYSAlPiUgCiAgICBmaWx0ZXIoc2V0ID09ICJNIikgJT4lIAogICAgc2VsZWN0KHBsYXNtaWQsIHN5bWJvbCwgdmFyID0ge3sgdmFyaWFibGUgfX0pICU+JSAKICAgIG11dGF0ZShmZyA9IHN0cl9zdWIoc3ltYm9sLCByZWdpb24sIHJlZ2lvbikgJT4lIHRvdXBwZXIoKSwKICAgICAgICAgICBiZyA9IHN5bWJvbCAlPiUgdG91cHBlcigpKQogICMgcmVwbGFjZSB0aGUgZm9yZWdyb3VuZCByZWdpb24gd2l0aCBYIGZvciBncm91cGluZwogIHN0cl9zdWIodG1wJGJnLCByZWdpb24sIHJlZ2lvbikgPC0gIlgiCiAgIyByZW9yZ2FuaXplIHRoZSB0aWJibGUgZm9yIGVhc2llciBoYW5kbGluZywgb3B0aW9uYWwKICB0bXAgPC0gcmVsb2NhdGUodG1wLCBmZywgYmcsIC5iZWZvcmUgPSBzeW1ib2wpICU+JSBzZWxlY3QoLXN5bWJvbCkKICAjIHBpdm90IHRoZSBkYXRhIGludG8gYSB3aWRlIGZvcm1hdCBzdWNoIHRoYXQgZm9yIGVhY2ggYmFja2dyb3VuZCwgdGhlcmUKICAjIGFyZSB0d28gdmFsdWVzIGZvciB0aGUgdmFyaWFibGUgb2YgaW50ZXJlc3QsIG9uZSBmcm9tIHRoZSBjaGltZXJhIHdpdGggCiAgIyBDZ1BobzQncyB2ZXJzaW9uIGluIHRoZSBmb3JlZ3JvdW5kIGFuZCBhbm90aGVyIHdpdGggU2NQaG80J3MgdmVyc2lvbgogIHRtcCA8LSB0bXAgJT4lIAogICAgc2VsZWN0KHBsYXNtaWQsIGZnLCBiZywgdmFyKSAlPiUgCiAgICBwaXZvdF93aWRlcihpZF9jb2xzID0gYmcsIG5hbWVzX2Zyb20gPSAiZmciLCAKICAgICAgICAgICAgICAgIHZhbHVlc19mcm9tID0gYyhwbGFzbWlkLCB2YXIpKSAlPiUgCiAgICB1bml0ZShwbGFzbWlkLCBzdGFydHNfd2l0aCgicGxhc21pZCIpKSAlPiUKICAgIG11dGF0ZShsYWJlbCA9IHBhc3RlKGJnLCBwbGFzbWlkLCBzZXAgPSAiXG4iKSkKICByZXR1cm4odG1wKQp9CmBgYAoKYGBgICAgICAgICAgCnggPSA1CnAxIDwtIG15X3Bsb3RfcmVnaW9uX2VmZmVjdF9vbmV2YXIoeCwgIkFfUEhPMiIpCnAyIDwtIG15X3Bsb3RfcmVnaW9uX2VmZmVjdF9vbmV2YXIoeCwgIkFfcGhvMiIpCnN1YnBsb3QocDEsIHAyLCBtYXJnaW4gPSAwLjA1KSAlPiUgCiAgbGF5b3V0KHRpdGxlID0gcGFzdGUoIlJlZ2lvbiIsIHgsICJzd2FwIGVmZmVjdCBvbiBBX1BITzIgYW5kIEFfcGhvMiIsIHNlcCA9ICIgIiksCiAgICAgICAgIHhheGlzID0gbGlzdCh0aXRsZSA9IHBhc3RlMCgiUmVnaW9uICIsIHgsICIgZnJvbSBDZ1BobzQiKSksCiAgICAgICAgIHlheGlzID0gbGlzdCh0aXRsZSA9IHBhc3RlMCgiUmVnaW9uICIsIHgsICIgZnJvbSBTY1BobzQiKSkgKQpgYGAKCkhlcmUsIEknZCBsaWtlIHRvIHRha2Ugd2hhdCBJIGJ1aWxkIGFib3ZlIGFuZCBjcmVhdGUgYSBuZXcgdGliYmxlLCBpbiB3aGljaCBlYWNoIHJvdyBpcyBhIGRpZmZlcmVudCBiYWNrZ3JvdW5kIChtYWtldXAgb2YgdGhlIGNoaW1lcmEgZXhjZXB0IGZvciB0aGUgZm9jYWwgcmVnaW9uKS4gVGhlIHZhbHVlIGNvbHVtbnMgYXJlOgoKMS4gIGRBX1BITzIgPSBBX1BITzJfQ2cgLSBBX1BITzJfU2MKMi4gIGRBX3BobzIgPSBBX3BobzJfQ2cgLSBBX3BobzJfU2MKMy4gIEFfUEhPMl9TYyA9IEFfUEhPMl9TYwoKVGhlIGdvYWwgaXMgdG8gcGxvdCBkQV9QSE8yIGFuZCBkQV9waG8yIHNpZGUtYnktc2lkZSBmb3IgZWFjaCBiYWNrZ3JvdW5kLgoKYGBge3J9Cm15X2NvbXBfcmVnaW9uX2VmZmVjdCA8LSBmdW5jdGlvbihyZWdpb24pewogICMgdGhpcyBmdW5jdGlvbiB1c2VzIG15X2NhbGNfcmVnaW9uX2VmZmVjdCB0byBnZXQgdGhlIHZhbHVlIGZvciB0aGUgdmFyaWFibGUgb2YgaW50ZXJlc3QKICAjIHdpdGggZWl0aGVyIENnIG9yIFNjIHZlcnNpb24gaW4gdGhlIGZvY2FsIHJlZ2lvbiwgc2VwYXJhdGVseSBmb3IgZWFjaCBiYWNrZ3JvdW5kIGNvbXBvc2l0aW9uCiAgIyBpdCBkb2VzIHNvIGZvciB0d28gdmFyaWFibGVzLCBBX1BITzIgYW5kIEFfcGhvMiwgdGhlbiBjYWxjdWxhdGUgZEFfUEhPMiwgZEFfcGhvMiwgYW5kCiAgIyBjb21iaW5lIHRoZW0KICBQSE8yID0gbXlfY2FsY19yZWdpb25fZWZmZWN0KHJlZ2lvbiwgIkFfUEhPMiIpICU+JSAKICAgIG11dGF0ZShkQV9QSE8yID0gdmFyX0MgLSB2YXJfUywKICAgICAgICAgICAjIG1lYW4gQV9QSE8yCiAgICAgICAgICAgTV9QSE8yID0gKHZhcl9TICsgdmFyX0MpLzIsCiAgICAgICAgICAgTkYgPSBpZmVsc2UoTV9QSE8yIDw9My41LCBUUlVFLCBGQUxTRSkpICU+JSAKICAgIHNlbGVjdCgtdmFyX1MsIC12YXJfQykKICAKICBwaG8yID0gbXlfY2FsY19yZWdpb25fZWZmZWN0KHJlZ2lvbiwgIkFfcGhvMiIpICU+JSAKICAgIG11dGF0ZShkQV9waG8yID0gdmFyX0MgLSB2YXJfUywgCiAgICAgICAgICAgTV9waG8yID0gKHZhcl9TICsgdmFyX0MpLzIpICU+JSAKICAgIHNlbGVjdCgtdmFyX1MsIC12YXJfQykKICAKICBkYXQgPC0gZnVsbF9qb2luKFBITzIsIHBobzIsIGJ5ID0gYygiYmciLCAicGxhc21pZCIsICJsYWJlbCIpKSAlPiUgCiAgICBzZWxlY3QoYmcsIHBsYXNtaWQsIGRBX1BITzIsIGRBX3BobzIsIE1fUEhPMiwgTV9waG8yLCBORikKICAKICByZXR1cm4oZGF0KQp9CgoKYGBgCgpgYGB7cn0KbXlfcGxvdF9yZWdpb25fZWZmZWN0X3R3b3Zhcl9saW5lKCIxIiwgIjQiKSMgJT4lIGdncGxvdGx5KCkKZ2dzYXZlKCIuLi9pbWcvMjAyNDAzMTAtcmVnaW9uLXN3YXAtZWZmZWN0LTEtb24tNC5wbmciLCB3aWR0aCA9IDYsIGhlaWdodCA9IDQpCm15X3Bsb3RfcmVnaW9uX2VmZmVjdF90d292YXJfbGluZSgiMyIsICI0IikjICU+JSBnZ3Bsb3RseSgpCmdnc2F2ZSgiLi4vaW1nLzIwMjQwMzEwLXJlZ2lvbi1zd2FwLWVmZmVjdC0zLW9uLTQucG5nIiwgd2lkdGggPSA2LCBoZWlnaHQgPSA0KQoKYGBgCgpgYGAgICAgICAgICAKbXlfcGxvdF9yZWdpb25fZWZmZWN0X3R3b3Zhcl9saW5lKCI0IiwgIjUiKSMgJT4lIGdncGxvdGx5KCkKZ2dzYXZlKCIuLi9pbWcvMjAyMzEyMjEtcmVnaW9uLXN3YXAtZWZmZWN0LTQtb24tNS5wbmciLCB3aWR0aCA9IDYsIGhlaWdodCA9IDQsIGRwaSA9IDE1MCkKbXlfcGxvdF9yZWdpb25fZWZmZWN0X3R3b3Zhcl9saW5lKCI1IiwgIjQiKSMgJT4lIGdncGxvdGx5KCkKZ2dzYXZlKCIuLi9pbWcvMjAyMzEyMjQtcmVnaW9uLXN3YXAtZWZmZWN0LTUtb24tNC5wbmciLCB3aWR0aCA9IDYsIGhlaWdodCA9IDQsIGRwaSA9IDIwMCkKYGBgCgpUaGUgbWFpbiBwbG90dGluZyBmdW5jdGlvbnMgYXJlIG5vdyBpbiBhIHNlcGFyYXRlIHNjcmlwdCBmaWxlIGluIGAuLi9zY3JpcHRgLiBUaGUgcGxvdHRpbmcgZnVuY3Rpb24gYmVsb3cgaXMgdG8gYWRhcHQgdGhlIHBsb3QgZm9yIGEgZmlndXJlIGluIHRoZSBwYXBlciwgc2ltdWx0YW5lb3VzbHkgc2hvd2luZyByZWdpb25zIDEtMy4KCmBgYHtyfQpteV9wbG90X3JlZ2lvbl9lZmZlY3RfdHdvdmFyX2xpbmVfcGFyIDwtIGZ1bmN0aW9uKHJlZ2lvbnMpewogICMgdGhpcyBmdW5jdGlvbiB1c2VzIG15X2NvbXBfcmVnaW9uX2VmZmVjdCB0byBnZW5lcmF0ZSB0aGUgZGF0YQogICMgYW5kIHBsb3QgdGhlIGRpZmZlcmVuY2UgaW4gQV9QSE8yIGFuZCBBX3BobzIgYmV0d2VlbiB0aGUgQ2dQaG80IHZzIFNjUGhvNAogICMgaW4gdGhlIGZvY2FsIHJlZ2lvbgogIGRhdCA8LSBtYXBfZGZyKHJlZ2lvbnMsIFwocmVnaW9uKSBteV9jb21wX3JlZ2lvbl9lZmZlY3QocmVnaW9uKSwgLmlkID0gInJlZ2lvbiIpICU+JSAKICAgIHBpdm90X2xvbmdlcihjb2xzID0gYyhkQV9QSE8yLCBkQV9waG8yKSwgCiAgICAgICAgICAgICAgICAgbmFtZXNfdG8gPSAiaG9zdCIsIHZhbHVlc190byA9ICJkaWZmIikgJT4lIAogICAgbXV0YXRlKGhvc3QgPSBmY3RfcmVjb2RlKGhvc3QsIGBQSE8yYCA9ICJkQV9QSE8yIiwgYHBobzLiiIZgID0gImRBX3BobzIiKSwKICAgICAgICAgICBob3N0ID0gZmN0X3JlbGV2ZWwoaG9zdCwgIlBITzIiKSkKICAjIHNwZWNpZnkgZ3JvdXBpbmcgdmFyaWFibGUKICBkYXQgPC0gbXV0YXRlKGRhdCwgCiAgICAgICAgICAgICAgICBncnAgPSBzdHJfc3ViKGJnLCA0LCA0KSAlPiUgdG91cHBlcigpLAogICAgICAgICAgICAgICAgZ3JwID0gZmN0X3JlY29kZShncnAsIENnUGhvNCA9ICJDIiwgU2NQaG80ID0gIlMiKSkjLAogICAgICAgICAgICAgICAgI3NoID0gc3RyX3N1YihiZywgNSwgNSkgJT4lIHRvdXBwZXIoKSwKICAgICAgICAgICAgICAgICNzaCA9IGZjdF9yZWNvZGUoc2gsIENnUGhvNCA9ICJDIiwgU2NQaG80ID0gIlMiKSApCiAgIyBzcGVjaWZ5IGFycm93IGFubm90YXRpb24KICBhcnJvdy54ID0gMC43CiAgYXJyb3cueSA9IChtYXgoZGF0JGRpZmYpIC0gbWluKGRhdCRkaWZmKSkgLyA1IAogICMgcGxvdAogIHAgPC0gZGF0ICU+JSAKICAgIGdncGxvdChhZXMoeCA9IGhvc3QsIHkgPSBkaWZmLCBsYWJlbCA9IGJnKSkgKwogICAgZ2VvbV9wb2ludChhZXMoY29sb3IgPSBncnApLCBzaXplID0gMiwgYWxwaGEgPSAwLjgsCiAgICAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25faml0dGVyKDAuMSkpICsgCiAgICBnZW9tX2xpbmUoYWVzKGdyb3VwID0gYmcpLCBsaW5ld2lkdGggPSAwLjIsIGFscGhhID0gMC44KSArCiAgICBmYWNldF9ncmlkKHJlZ2lvbiB+IGdycCwgbGFiZWxsZXIgPSBsYWJlbGxlcigKICAgICAgZ3JwID0gYyhDZ1BobzQgPSAiUDJJRDpDZyIsIFNjUGhvNCA9ICJQMklEOlNjIiksCiAgICAgIHJlZ2lvbiA9IGxhYmVsX2JvdGgKICAgICkpICsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCgiUDJJRDoiLCB2YWx1ZXMgPSBjKCJvcmFuZ2UiLCAiZ3JheTMwIiksIGd1aWRlID0gIm5vbmUiKSArCiAgICAjc2NhbGVfc2hhcGVfbWFudWFsKCJEQkQ6IiwgdmFsdWVzID0gYygxOSwgMSkpICsKICAgIHlsYWIoIlJlZ2lvbiBzd2FwIGVmZmVjdCAoQ2ctU2MpIikgKwogICAgdGhlbWVfYncoYmFzZV9zaXplID0gMTgpICsgCiAgICB0aGVtZSgKICAgICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF9ibGFuaygpLAogICAgICBheGlzLnRpdGxlLnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgwLjkpKSwKICAgICAgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoZmFjZSA9IDMpLAogICAgICBheGlzLnRleHQueSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDAuOCkpLAogICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDAuOCkpLAogICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgwLjkpKSwKICAgICAgbGVnZW5kLnBvc2l0aW9uID0gInRvcCIsCiAgICAgIHN0cmlwLmJhY2tncm91bmQgPSBlbGVtZW50X2JsYW5rKCkKICAgICkKICByZXR1cm4ocCkKfQpgYGAKCmBgYHtyfQpteV9wbG90X3JlZ2lvbl9lZmZlY3RfdHdvdmFyX2xpbmVfcGFyKGMoMSwzKSkKZ2dzYXZlKCIuLi9pbWcvMjAyNDAzMTAtcmVnaW9uLXN3YXAtZWZmZWN0LTFuMy1vbi00LnBuZyIsCiAgICAgICB3aWR0aCA9IDUsIGhlaWdodCA9IDMuNSkKYGBgCgpgYGB7cn0KbXlfcGxvdF9yZWdpb25fZWZmZWN0X3R3b3Zhcl9zaWRlIDwtIGZ1bmN0aW9uKHJlZ2lvbil7CiAgIyB0aGlzIGZ1bmN0aW9uIHVzZXMgbXlfY29tcF9yZWdpb25fZWZmZWN0IHRvIGdlbmVyYXRlIHRoZSBkYXRhCiAgIyBhbmQgcGxvdCB0aGUgZGlmZmVyZW5jZSBpbiBBX1BITzIgYW5kIEFfcGhvMiBiZXR3ZWVuIHRoZSBDZ1BobzQgdnMgU2NQaG80CiAgIyBpbiB0aGUgZm9jYWwgcmVnaW9uCiAgZGF0IDwtIG15X2NvbXBfcmVnaW9uX2VmZmVjdChyZWdpb24pICU+JSAKICAgIHBpdm90X2xvbmdlcihjb2xzID0gYyhkQV9QSE8yLCBkQV9waG8yKSwgCiAgICAgICAgICAgICAgICAgbmFtZXNfdG8gPSAiaG9zdCIsIHZhbHVlc190byA9ICJkaWZmIikgJT4lIAogICAgbXV0YXRlKGhvc3QgPSBmY3RfcmVjb2RlKGhvc3QsIGBQSE8yYCA9ICJkQV9QSE8yIiwgYHBobzJgID0gImRBX3BobzIiKSwKICAgICAgICAgICBob3N0ID0gZmN0X3JlbGV2ZWwoaG9zdCwgIlBITzIiKSkKICBwIDwtIGRhdCAlPiUgCiAgICBnZ3Bsb3QoYWVzKHggPSBiZywgeSA9IGRpZmYsIGdyb3VwID0gaG9zdCkpICsKICAgIGdlb21fY29sKGFlcyhmaWxsID0gaG9zdCksIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2UoMC45KSkgKwogICAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gaG9zdC5jb2xvcnMpICsKICAgIHlsYWIoIlJlZ2lvbiBzd2FwIGRpZmYgKENnIHZzIFNjKSIpICsKICAgIHRoZW1lX2Nvd3Bsb3QoZm9udF9zaXplID0gMjApICsgCiAgICBwYW5lbF9ib3JkZXIoY29sb3IgPSAiZ3JheTMwIikgKwogICAgYmFja2dyb3VuZF9ncmlkKG1ham9yID0gInkiLCBtaW5vciA9ICJub25lIikgKwogICAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgaGp1c3QgPSAxLCBmYW1pbHkgPSAiY291cmllciIpLAogICAgICAgICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgICAgbGVnZW5kLnBvc2l0aW9uID0gInRvcCIpCiAgcmV0dXJuKHApCn0KYGBgCgojIyBQMklEOkNnX0RCRDpTYyBmYWlsCgpIaWdobGlnaHQgdGhlIHN1YnNldCBvZiB0aGUgY2hpbWVyYXMgd2l0aCBQMklEOkNnICsgREJEOlNjLCBtb3N0IG9mIHdoaWNoIGFyZSBub24gZnVuY3Rpb25hbAoKYGBge3J9Cm15X3NjYXR0ZXJfcGxvdCgiWFhYWENTIikgKwogICAgbGFicyh4ID0gYnF1b3RlKEFbUEhPMl0pLCB5ID0gYnF1b3RlKEFbcGhvMl0pKQpnZ3NhdmUoZmlsZW5hbWUgPSAiLi4vaW1nLzIwMjQwMjE0LVBobzQtY2hpbWVyYXMtc2NhdHRlci1QMklEX0NnLURCRF9TYy5wbmciLAogICAgICAgcGxvdCA9IHAsIHdpZHRoID0gNiwgaGVpZ2h0ID0gNCwgZHBpID0gMzAwKQpgYGAKCmBgYHtyfQp4IDwtIG15X2RhdGFfc2VsZWN0KHBhdHRlcm4gPSAiWFhYWENTIiwgU2V0ID0gIk0iKQpteV9kYXRhX3ByZXAoeCkgJT4lIAogIG11dGF0ZShncm91cCA9IGZjdF9yZWNvZGUoZ3JvdXAsICJjaGltZXJhIiA9ICJuLmYuIikpICU+JSAKICBteV9wbG90X2NvbXBvbmVudHMoKQpnZ3NhdmUoIi4uL2ltZy8yMDI0MDIxMy1QMklEX0NnLURCRF9TYy1jb21wb25lbnRzLnBuZyIsIHdpZHRoID0gOCwgaGVpZ2h0ID0gNSkKYGBgCgojIyBUcmlhbmdsZSBoZWF0bWFwCgpGaXJzdCwgd3JpdGUgYSBmdW5jdGlvbiB0byBnZW5lcmF0ZSB0aGUgZGF0YSBmb3IgcGxvdHRpbmcuIElmIHdlIGFyZSBnb2luZyB0byB1c2UgZ2dwbG90LCB3ZSBuZWVkIGEgdGliYmxlIHRvIHN0b3JlIHRoZSBkYXRhLCBzb21ldGhpbmcgaW4gdGhlIGZvbGxvd2luZyBmb3JtCgp8IHBsYXNtaWQgfCBzeW1ib2wgfCBSZWdpb25BIHwgUmVnaW9uQiB8IEFfUEhPMiB8IEFfcGhvMiB8IHJBX1BITzIgfCBib29zdCB8IHBlcmNfcGhvMiB8Cnw6LS0tLS0tLS18Oi0tLS0tLS18Oi0tLS0tLS0tfDotLS0tLS0tLXw6LS0tLS0tLXw6LS0tLS0tLXw6LS0tLS0tLS18Oi0tLS0tLXw6LS0tLS0tLS0tLXwKfCAyMDkgICAgIHwgQ0NTQ0MgIHwgMyAgICAgICB8IDMgICAgICAgfCA4LjI1ICAgfCA3LjgyICAgfCAwLjQ2OCAgIHwgMS4wNiAgfCAwLjk0ICAgICAgfAoKSWYgd2UgYXJlIG9rIHdpdGggdXNpbmcgbm9uIGdncGxvdCAtIGhlYXRtYXBzIGFyZSBub3QgZ2dwbG90J3Mgc3RyZW5ndGggYW55d2F5cyAtIHdlIGNhbiBqdXN0IGJ1aWxkIGEgbWF0cml4LgoKTm90ZSB0aGF0IHRoaXMgd2F5IG9mIHN1bW1hcml6aW5nIHRoZSBkYXRhIGhhcyBtYW55IGxpbWl0YWl0b25zOiAxKSBpdCByZXF1aXJlcyBzcGVjaWZ5aW5nIHRoZSByZWZlcmVuY2UsIGVpdGhlciBDQ0NDQyBvciBTU1NTUy4gRXZlcnl0aGluZyBpcyBtZWFzdXJlZCBhZ2FpbnN0IHRoYXQ7IDIpIGl0IG9ubHkgc2hvd3MgcGFpcndpc2UgKHR3byByZWdpb24pIGludGVyYWN0aW9ucy4gVGhpcyB0dXJucyBvdXQgdG8gYmUgZmluZSB3aXRoIGZpdmUgcmVnaW9ucywgc2luY2UgZXZlcnkgY2hpbWVyYSBjYW4gYmUgZXhwcmVzc2VkIGFzIGVpdGhlciBhIDAsIDEgb3IgMiByZWdpb24gc3dhcCBmcm9tIG9uZSBvZiB0aGUgdHdvIHJlZmVyZW5jZSBnZW5vdHlwZXMuIFdpdGggNiBvciBtb3JlIHJlZ2lvbnMsIGhpZ2hlciBsZXZlbCAoMyBvciBtb3JlIHJlZ2lvbikgaW50ZXJhY3Rpb25zIGNhbm5vdCBiZSB2aXN1YWxpemVkIHRoaXMgd2F5LiBCZWNhdXNlIG9mIHRoaXMsIHdlIHdpbGwgZm9jdXMgb24ganVzdCB0aGUgbWFpbiBzZXQgZm9yIHRoaXMgYW5hbHlzaXMuCgpUbyBidWlsZCB0aGUgbWF0cml4LCB3ZSBuZWVkIHRvIGZpcnN0IGlkZW50aWZ5IHRoZSBjaGltZXJhcyB0aGF0IGJlbG9uZyB0byB0aGUgc2V0LiBGb3IgdGhhdCwgd2Ugd2lsbCB1c2UgdGhlICJtYWluIiBzZXQsIHdpdGggdGhlIGZpdmUgcmVnaW9uIHNwbGl0LCBmb3IgdGhlIG1vbWVudCBhdCBsZWFzdC4gVGhlIGZ1bmN0aW9uIHdpbGwgZmlyc3QgZGV0ZXJtaW5lIHdoaWNoIHJlZmVyZW5jZSB0byB1c2UuIElmIHdlIHVzZSBTU1NTUyBhcyB0aGUgcmVmZXJlbmNlLCBmb3IgZXhhbXBsZSwgd2Ugd2lsbCBhc3NpZ24gMCB0byB0aGUgcmVmZXJlbmNlLiBBbGwgb3RoZXIgY2hpbWVyYXMgd2l0aCAxIG9yIDIgcmVnaW9ucyBmcm9tIENnIHdpbGwgYmUgdXNlZCB0byBmaWxsIGFuIHVwcGVyIHRyaWFuZ3VsYXIgbWF0cml4LCB1c2luZyBvbmUgb2YgdGhlIHZhbHVlcyBvZiBpbnRlcmVzdCwgZS5nLiwgQV9QSE8yLgoKYGBge3J9Cm15X3VwcGVyX3RyaWFuZ3VsYXJfbWF0IDwtIGZ1bmN0aW9uKGFsdCA9ICJDIiwgdmFyID0gIkFfUEhPMiIsIG5mLmFzLm5hID0gRil7CiAgIyBnaXZlbiB0aGUgYWx0ZXJuYXRpdmUgYWxsZWxlIChDL1MpIGFuZCBhIHZhcmlhYmxlIG9mIGludGVyZXN0LCBlLmcuLCBBX1BITzIsCiAgIyBvdXRwdXQgYW4gdXBwZXIgdHJpYW5ndWxhciBtYXRyaXggY29udGFpbmluZyB0aGUgdmFsdWVzIGZyb20gdGhlIHZhcmlhYmxlIAogICMgb2YgaW50ZXJlc3QsIHdpdGggdGhlIHJvdyBhbmQgY29sIG51bWJlcnMgYmFzZWQgb24gdGhlIGZpcnN0IGFuZCBzZWNvbmQKICAjIHBvc2l0aW9ucyBjb250YWluaW5nIHRoZSBhbHRlcm5hdGl2ZSBhbGxlbGUuIElmIGFsbCBwb3NpdGlvbnMgY29udGFpbiB0aGUgCiAgIyByZWZlcmVuY2UgYWxsZWxlLCB0aGUgdmFsdWUgaXMgc3VidHJhY3RlZCBmcm9tIGFsbCB2YWx1ZXMgaW4gdGhlIG1hdHJpeAogICMgd2hlbiBqdXN0IG9uZSBwb3NpdGlvbiBpcyB0aGUgYWx0ZXJuYXRpdmUgYWxsZWxlLCB0aGUgdmFsdWUgaW4gdGhlIGRpYWdvbmFsCiAgIyBpcyBzZXQuIHdoZW4gdGhlcmUgYXJlIG1vcmUgdGhhbiAyIHJlZ2lvbnMgY29udGFpbmluZyB0aGUgYWx0ZXJuYXRpdmUgYWxsZWxlCiAgIyBza2lwLgogICMgaWYgIm5mLmFzLm5hID0gVFJVRSIsIGV2YWx1YXRlIGlmIHRoZSBhY3Rpdml0eSBvZiBlaXRoZXIgb2YgdGhlIHR3byBjaGltZXJhcwogICMgYmVpbmcgY29tcGFyZWQgaXMgbm9uIGZ1bmN0aW9uYWwuIGlmIHllcywgc2V0IHRoZSBjb3JyZXNwb25kaW5nIG1hdHJpeCB2YWx1ZQogICMgdG8gTkEKICBvdXRfbWF0IDwtIG1hdHJpeChOQSwgbnJvdyA9IDUsIG5jb2wgPSA1KQogIHJlZl92YWwgPC0gTkEKICBkYXQgPC0gZmlsdGVyKHhpbWVyYSwgc2V0ID09ICJNIikgJT4lIAogICAgbXV0YXRlKFMgPSBhcy5jaGFyYWN0ZXIoc3ltYm9sKSAlPiUgdG91cHBlcigpKQogIGlmKG5mLmFzLm5hKXsKICAgIGRhdCA8LSBmaWx0ZXIoZGF0LCBncm91cCAhPSAibi5mLiIpCiAgfQogIGZvcihpIGluIHNlcSgxLCBucm93KGRhdCkpKXsKICAgIHN5bWJvbCA9IGRhdFtpLCAiUyJdCiAgICAjIGRldGVybWluZSB3aGljaCBwb3NpdGlvbnMgY29udGFpbiB0aGUgYWx0ZXJuYXRpdmUgYWxsZWxlCiAgICBwID0gc3RyX2xvY2F0ZV9hbGwoc3ltYm9sLCBhbHQpW1sxXV1bLCJzdGFydCJdCiAgICBsID0gbGVuZ3RoKHApICAgIyBob3cgbWFueSBwb3NpdGlvbnMgY29udGFpbiB0aGUgYWx0IGFsbGVsZQogICAgdiA9IGRhdFtbdmFyXV1baV0gIyByZXRyaWV2ZSB0aGUgdmFsdWUgb2YgdGhlIHZhcmlhYmxlCiAgICBpZihsID09IDApCiAgICAgIHJlZl92YWwgPSB2CiAgICBlbHNlIGlmKGwgPT0gMSkKICAgICAgb3V0X21hdFtwLCBwXSA9IHYKICAgIGVsc2UgaWYobCA9PSAyKQogICAgICBvdXRfbWF0W3BbMV0sIHBbMl1dID0gdgogIH0KICBvdXRfbWF0ID0gb3V0X21hdCAtIHJlZl92YWwKICByZXR1cm4ob3V0X21hdCkKfQpgYGAKCmBgYHtyfQpteV9jb21iaW5lZF90cmlhbmd1bGFyX21hdCA8LSBmdW5jdGlvbihhbHQgPSAiQyIpewogICMgZ2l2ZW4gdGhlIGFsdGVybmF0aXZlIGFsbGVsZSAoQy9TKSwgb3V0cHV0IGEgbWF0cml4IGNvbnRhaW5pbmcgdGhlIHZhbHVlcwogICMgZm9yIGJvdGggd2l0aCBhbmQgd2l0aG91dCBQaG8yLCBhcnJhbmdlZCBpbiB0d28gY29tcGxlbWVudGFyeSB0cmlhZ3VsYXIKICAjIG1hdHJpY2VzLCB3aXRoIHRoZSByb3cgYW5kIGNvbCBudW1iZXJzIGJhc2VkIG9uIHRoZSBmaXJzdCBhbmQgc2Vjb25kCiAgIyBwb3NpdGlvbnMgY29udGFpbmluZyB0aGUgYWx0ZXJuYXRpdmUgYWxsZWxlLiBJZiBhbGwgcG9zaXRpb25zIGNvbnRhaW4gdGhlIAogICMgcmVmZXJlbmNlIGFsbGVsZSwgdGhlIHZhbHVlIGlzIHN1YnRyYWN0ZWQgZnJvbSBhbGwgdmFsdWVzIGluIHRoZSBtYXRyaXgKICAjIHdoZW4ganVzdCBvbmUgcG9zaXRpb24gaXMgdGhlIGFsdGVybmF0aXZlIGFsbGVsZSwgdGhlIHZhbHVlIGluIHRoZSBkaWFnb25hbAogICMgaXMgc2V0LiB3aGVuIHRoZXJlIGFyZSBtb3JlIHRoYW4gMiByZWdpb25zIGNvbnRhaW5pbmcgdGhlIGFsdGVybmF0aXZlIGFsbGVsZQogICMgc2tpcC4KICBvdXRfbWF0IDwtIG1hdHJpeChOQSwgbnJvdyA9IDYsIG5jb2wgPSA2KQogIHVwcGVyIDwtIGNiaW5kKE5BLCBteV91cHBlcl90cmlhbmd1bGFyX21hdChhbHQsIHZhciA9ICJBX1BITzIiLCApKSAlPiUgCiAgICByYmluZCguLCBOQSkKICBsb3dlciA8LSByYmluZChOQSwgdChteV91cHBlcl90cmlhbmd1bGFyX21hdChhbHQsIHZhciA9ICJBX3BobzIiKSkpICU+JSAKICAgIGNiaW5kKC4sIE5BKQogIG91dF9tYXQgPSBpZmVsc2UoaXMubmEodXBwZXIpLCBsb3dlciwgdXBwZXIpCiAgcmV0dXJuKG91dF9tYXQpCn0KYGBgCgpgYGB7cn0KbXlfcGxvdF90cmlhbmdsZV9oZWF0bWFwIDwtIGZ1bmN0aW9uKGFsdCwgdmFyKXsKICAjIHRoaXMgZnVuY3Rpb24gdGFrZXMgdGhlIG91dHB1dCBvZiB0aGUgZnVuY3Rpb24gYWJvdmUgYW5kIG1ha2VzIGEgaGVhdG1hcAogICMgdXNpbmcgcGhlYXRtYXAgZnVuY3Rpb24sIHRoZW4gcm90YXRlcyBpdCB1c2luZyBncmlkIGdyYXBoaWNzCiAgIyB0aGFua3MgdG8gaHR0cHM6Ly9ib29rZG93bi5vcmcvcmRwZW5nL1JQcm9nREEvdGhlLWdyaWQtcGFja2FnZS5odG1sI2dyaWQtZ3JhcGhpY3MtY29vcmRpbmF0ZS1zeXN0ZW1zCiAgIyBhZGRpbmcgdGl0bGUgYmFzZWQgb24gaHR0cHM6Ly9kYXZldGFuZy5naXRodWIuaW8vbXVzZS9waGVhdG1hcC5odG1sCiAgCiAgIyBjb25zdHJ1Y3QgdGl0bGUgb2YgcGxvdAogIHJlZiA9IGlmZWxzZShhbHQgPT0gIkMiLCAiU2NQaG80IiwgIkNnUGhvNCIpCiAgYmcgPSBpZmVsc2UodmFyID09ICJBX1BITzIiLCAid2l0aCBQSE8yIiwgIncvbyBwaG8yIikKICBteV90aXRsZSA8LSBwYXN0ZSgiRXBpc3Rhc2lzIGJldHdlZW4gcmVnaW9ucyBvbiIsIHJlZiwgImJhY2tncm91bmQiLCBiZykKICB0ZXN0IDwtIG15X3VwcGVyX3RyaWFuZ3VsYXJfbWF0KGFsdCA9IGFsdCwgdmFyID0gdmFyKQogIHBhbGV0dGVMZW5ndGggPSA1MAogIG15Q29sb3JzIDwtIGNvbG9yUmFtcFBhbGV0dGUoYygic3RlZWxibHVlMyIsICJncmF5OTAiLCAicmVkIikpKHBhbGV0dGVMZW5ndGgpCiAgcm5nIDwtIG1heChhYnModGVzdCksIG5hLnJtID0gVFJVRSkKICBteUJyZWFrcyA8LSBjKHNlcSgtcm5nLCAwLCBsZW5ndGgub3V0PWNlaWxpbmcocGFsZXR0ZUxlbmd0aC8yKSArIDEpLCAKICAgICAgICAgICAgICAgIHNlcShybmcvcGFsZXR0ZUxlbmd0aCwgcm5nLAogICAgICAgICAgICAgICAgICAgIGxlbmd0aC5vdXQ9Zmxvb3IocGFsZXR0ZUxlbmd0aC8yKSkpCiAgcCA8LSBwaGVhdG1hcDo6cGhlYXRtYXAodGVzdCwgY29sb3IgPSBteUNvbG9ycywgYnJlYWtzID0gbXlCcmVha3MsCiAgICAgICAgICAgICAgICAgICAgICAgICAgYm9yZGVyX2NvbG9yID0gTkEsIG5hX2NvbCA9IE5BLCBzaWxlbnQgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICAgICAgIGNsdXN0ZXJfY29scyA9IEZBTFNFLCBjbHVzdGVyX3Jvd3MgPSBGQUxTRSkKICB2cCA8LSB2aWV3cG9ydCh4ID0gMC41LCB5ID0gMC4yNSwKICAgICAgICAgICAgICAgICB3aWR0aCA9IHVuaXQoNC41LCAiaW4iKSwgaGVpZ2h0ID0gdW5pdCg0LjUsICJpbiIpLCBhbmdsZSA9IDQ3KSAKICBncmlkLm5ld3BhZ2UoKQogIHB1c2hWaWV3cG9ydCh2cCkKICBncmlkLmRyYXcocCRndGFibGUpCiAgcG9wVmlld3BvcnQoKQogIGdyaWQudGV4dChsYWJlbCA9IG15X3RpdGxlLCB4ID0gMC41LCB5ID0gMC45NSwgZ3AgPSBncGFyKGZvbnRzaXplID0gMTYsIGZvbnRmYWNlID0gImJvbGQiKSkKICByZXR1cm4ocCkKfQpgYGAKCmBgYHtyfQpteV9wbG90X2NvbWJpbmVkX3RyaWFuZ2xlX2hlYXRtYXAgPC0gZnVuY3Rpb24oYWx0KXsKICAjIHRoaXMgZnVuY3Rpb24gdGFrZXMgdGhlIG91dHB1dCBvZiB0aGUgZnVuY3Rpb24gbXlfY29tYmluZWRfdHJpYW5ndWxhcl9tYXQoKQogICMgdXNpbmcgcGhlYXRtYXAgZnVuY3Rpb24sIHRoZW4gcm90YXRlcyBpdCB1c2luZyBncmlkIGdyYXBoaWNzCiAgIyB0aGFua3MgdG8gaHR0cHM6Ly9ib29rZG93bi5vcmcvcmRwZW5nL1JQcm9nREEvdGhlLWdyaWQtcGFja2FnZS5odG1sI2dyaWQtZ3JhcGhpY3MtY29vcmRpbmF0ZS1zeXN0ZW1zCiAgIyBhZGRpbmcgdGl0bGUgYmFzZWQgb24gaHR0cHM6Ly9kYXZldGFuZy5naXRodWIuaW8vbXVzZS9waGVhdG1hcC5odG1sCiAgCiAgIyBjb25zdHJ1Y3QgdGl0bGUgb2YgcGxvdAogIHJlZiA9IGlmZWxzZShhbHQgPT0gIkMiLCAiU2NQaG80IiwgIkNnUGhvNCIpCiAgbXlfdGl0bGUgPC0gcGFzdGUoIkVwaXN0YXNpcyBiZXR3ZWVuIHJlZ2lvbnMgb24iLCByZWYsICJiYWNrZ3JvdW5kIikKICB0ZXN0IDwtIG15X2NvbWJpbmVkX3RyaWFuZ3VsYXJfbWF0KGFsdCA9IGFsdCkKICBwYWxldHRlTGVuZ3RoID0gNTAKICBteUNvbG9ycyA8LSBjb2xvclJhbXBQYWxldHRlKGMoInN0ZWVsYmx1ZSIsICJncmF5OTAiLCAicmVkIikpKHBhbGV0dGVMZW5ndGgpCiAgcm5nIDwtIG1heChhYnModGVzdCksIG5hLnJtID0gVFJVRSkKICBteUJyZWFrcyA8LSBjKHNlcSgtcm5nLCAwLCBsZW5ndGgub3V0PWNlaWxpbmcocGFsZXR0ZUxlbmd0aC8yKSArIDEpLCAKICAgICAgICAgICAgICAgIHNlcShybmcvcGFsZXR0ZUxlbmd0aCwgcm5nLAogICAgICAgICAgICAgICAgICAgIGxlbmd0aC5vdXQ9Zmxvb3IocGFsZXR0ZUxlbmd0aC8yKSkpCiAgcCA8LSBwaGVhdG1hcDo6cGhlYXRtYXAodGVzdCwgY29sb3IgPSBteUNvbG9ycywgYnJlYWtzID0gbXlCcmVha3MsCiAgICAgICAgICAgICAgICAgICAgICAgICAgYm9yZGVyX2NvbG9yID0gTkEsIG5hX2NvbCA9IE5BLCBzaWxlbnQgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICAgICAgIGNsdXN0ZXJfY29scyA9IEZBTFNFLCBjbHVzdGVyX3Jvd3MgPSBGQUxTRSkKICB2cCA8LSB2aWV3cG9ydCh4ID0gMC41LCB5ID0gMC40NSwKICAgICAgICAgICAgICAgICB3aWR0aCA9IHVuaXQoMywgImluIiksIGhlaWdodCA9IHVuaXQoMi44LCAiaW4iKSwgYW5nbGUgPSA0NykgCiAgZ3JpZC5uZXdwYWdlKCkKICBwdXNoVmlld3BvcnQodnApCiAgZ3JpZC5kcmF3KHAkZ3RhYmxlKQogIHBvcFZpZXdwb3J0KCkKICBncmlkLnRleHQobGFiZWwgPSBteV90aXRsZSwgeCA9IDAuNSwgeSA9IDAuOTUsIAogICAgICAgICAgICBncCA9IGdwYXIoZm9udHNpemUgPSAxNiwgZm9udGZhY2UgPSAiYm9sZCIpKQogIGdyaWQudGV4dChsYWJlbCA9ICJXaXRoIFBobzIiLCB4ID0gMC4xLCB5ID0gMC42NSwganVzdCA9IGMoImxlZnQiLCAidG9wIiksCiAgICAgICAgICAgIGdwID0gZ3Bhcihmb250c2l6ZSA9IDE0LCBmb250ZmFjZSA9ICJib2xkIikpCiAgZ3JpZC50ZXh0KGxhYmVsID0gIldpdGhvdXQgcGhvMiIsIHggPSAwLjEsIHkgPSAwLjI1LCBqdXN0ID0gYygibGVmdCIsICJ0b3AiKSwgCiAgICAgICAgICAgIGdwID0gZ3Bhcihmb250c2l6ZSA9IDE0LCBmb250ZmFjZSA9ICJib2xkIikpCiAgcmV0dXJuKHApCn0KYGBgCgpgYGB7cn0KcG5nKCIuLi9pbWcvMjAyNDAxMTUtdHJpYW5nbGUtaGVhdG1hcC1DZ1BobzQtcmVmLnBuZyIsIHdpZHRoID0gNywgaGVpZ2h0ID0gNSwgdW5pdHMgPSAiaW4iLCByZXMgPSAzMDApCnAxIDwtIG15X3Bsb3RfY29tYmluZWRfdHJpYW5nbGVfaGVhdG1hcCgiQyIpCmRldi5vZmYoKQpwbmcoIi4uL2ltZy8yMDI0MDExNS10cmlhbmdsZS1oZWF0bWFwLVNjUGhvNC1yZWYucG5nIiwgd2lkdGggPSA3LCBoZWlnaHQgPSA1LCB1bml0cyA9ICJpbiIsIHJlcyA9IDMwMCkKcDIgPC0gbXlfcGxvdF9jb21iaW5lZF90cmlhbmdsZV9oZWF0bWFwKCJTIikKZGV2Lm9mZigpCmBgYAo=